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The long-time asymptotics is analyzed for finite energy solutions of the ID Schrodinger 
equation coupled to a nonlinear oscillator; mathematically the system under study is a 

^^ ■ nonlinear Schrodinger equation, whose nonlinear term includes a Dirac delta. The coupled 

system is invariant with respect to the phase rotation group U{1). This article, which 
extends the results of a previous one, provides a proof of asymptotic stability of solitary 
wave solutions in the case that the linearization contains a single discrete oscillatory mode 

K> , satisfying a non-degeneracy assumption of the type known as the Fermi Golden Rule. 
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1 Introduction and statement of results 

1.1 Introduction 

In this article we continue the study of large time asymptotics for a model ?7(l)-invariant 
nonlinear Schrodinger equation 

iip(x,t)=-ip"{x,t)-5{x)F(ip(0,t)), xeR, (1.1) 

which was initiated in pp. Here ip(x,t) is a continuous complex- valued wave function and 
F is a continuous function, the dots stand for the derivatives in t and the primes in x; all 
derivatives and the equation are understood in the distribution sense. Our main focus is on the 
role that certain solitary waves (also referred to as nonlinear bound states, or solitons) play in 
the description of the solution for large times. These solitary waves are solutions of the form 
e' luJt ip U j(x) , where ip u solves a nonlinear eigenvalue problem ( 11.101) . In [1] the asymptotic stability 
of these solitary waves was proved under a condition on the nonlinearity which ensures that the 
linearization about the solitary wave consists entirely of continuous spectrum, except for the 
two dimensional generalized null space which is always present due to the U(l) symmetry of the 
equation. In this article this result is extended to the case that the spectrum of the linearization 
includes an additional discrete component, which satisfies a non- degeneracy condition related to 
the Fermi Golden Rule. In order to explain these results fully we will introduce our conditions 
on the nonlinearity F, in the remainder of this section. In the following section we will discuss 
the basic properties of the solitary waves. We will then be able to state our main theorem 
precisely as theorem 11.31 For a more lengthy discussion of our motivation, and of previous 
results in the literature ([2], [6], [TJ, [HI [El E]) we refer the reader to the introduction of [I]. 

It will be convenient to rewrite ( 11.11) in real form: we identify a complex number ip = ipi+iip2 
with the real two-dimensional vector (0i, ^2) £ R 2 and rewrite ( 11.11) in the vectorial form 

jj>{x, t) = -4"{x, t) - 5(x)F(i>(0, t)), (1.2) 

where 

';-;)■ (i - 3) 

and F(ip) G R 2 is the real vector version of F(ip) G C. We assume that the oscillator force F 
admits a real-valued potential 

F(V>) = -W(V>), fe« 2 , U G C 2 (R 2 ). (1.4) 

Then (II -2D is formally a Hamiltonian system with Hamiltonian 

n^) = l -fw\ 2 dx + u{^)) (1.5) 

which is conserved for sufficiently regular finite energy solutions. We assume that the potential 
U{ip) satisfies the inequality 

U(z)>A-B\z\ 2 with some AeR, B > 0. (1.6) 

We also assume that U(ip) = u(\ip\ 2 ) with u G C 2 (R). Therefore, by (Ti~4|) . 

F(V>) = a(\ij\ 2 )tp, ^GC, aeC^R), (1.7) 



where a(\ip\ 2 ) is real. Then F(e id ip) = e i9 F(ip), 9 G [0, 2vr], and e id ip(x, t) is a solution to (TO) 
if ijj(x,t) is. Therefore, equation (11.11) is [/(l)-invariant and the Nother theorem implies the 
conservation of the following charge: 

Q(iP) = / \ip\ 2 dx = const . (1.8) 

Under these conditions the existence of global solutions to the Cauchy problem for (11.11) 
was proved in [I]. We work in the space iJ 1 (R) = H 1 , the Sobolev space of complex valued 

measurable functions with / (|^'| 2 + |^| 2 )^ < oo, and C&(1R., X) is the space of bounded 

continuous functions M — > X into a Banach space X. 



Theorem 1.1 ([!]). Under conditions ( |1.^| ), (GUI' and ( [1.7| ) ; the following statements hold. 



i) For any ipo(x) G H 1 there exists a unique solution ip(t) = ip(- ,t) G Cb(R, H 1 ) to the 

equation U.l\) with initial condition ip(x,0) = 'tp^x). 

ii) The charge Q{ip{t)) and Hamiltonian T-C(tp(t)) are constant along the solution. 

Hi) There exists A(^o) > such that the following a priori bound holds: 

sup ||^(t) \\ H i <A(^ )<oo. (1.9) 

Next, in §1.21 we describe all nonzero solitary waves, and then formulate the main theorem 

mra 



1.2 Solitary waves 

Equation (11.11) admits finite energy solutions of type ip u} (x)e tuJt , called solitary waves or non- 
linear eig en functions. The frequency u and the amplitude ip^x) solve the following nonlinear 
eigenvalue problem: 

- uil> u (x) = -£(x) - 5(x)F(ij> u (0)), xeR. (1.10) 

It is straightforward to check (see [T]) that the set of all nonzero solitary waves consists of 
functions ip ul (x)e t9 , ipoj(x) = C{uj)e~^ x \ C > 0, u > 0, where 

yfc = a(C 2 )/2 > 0, 

and where 9 G [0, 2tt] is arbitrary. This condition means that C is restricted to lie in a set 
which, in the case of polynomial F, is a finite union of one-dimensional intervals. Notice that 
C = corresponds to the zero function ip(x) = which is always a solitary wave as F(0) = 0, 
and for uj < only the zero solitary wave exists. The real form of the solitary wave is e^^ 
where ^ = (^(x), 0). We will also need the following lemma from pQ: 

Lemma 1.2. For C > 0, a > we have 

du / |^(x)| 2 cfc > if < a' < a/C 2 . 



Linearization at the solitary wave e- ?6l \l/ w leads to the operator 

B = -^- 2 +u-5(xMC 2 )+2a'(C 2 )C 2 P l ]=^ B Q 1 ^ ), (1.11) 

where P\ is the projector in R 2 acting as I I i— > I 

d 2 d 2 

Di = -— + u - S(x)[a + 2a C 2 }, D 2 = -j^ + to - S(x)a. 

(see PQ). Let C = J -1 B. The continuous spectrum of C coincides with C + U C_ where C + = 
[iuj,ioo), and C_ = (— zoo, — iuj}. The discrete spectrum depend is described in [TJ. Zero is 
always present in the discrete spectrum on account of the circular symmetry of the problem, 
and there is a corresponding generalized eigenspace of dimension at least two. If a' > a/C 2 
there is a positive eigenvalue corresponding to linear instability of the solitary wave, while for 
a' < a/C 2 the discrete spectrum consists either only of zero, or contains in addition two pure 
imaginary eigenvalues. The presence of such imaginary discrete spectrum corresponds to the 
possibility of a periodic pulsation of the solitary wave at the linearized level, a possibility which 
has to be taken care of in the proof of asymptotic stability. 

1.3 Statement of the main theorem 

Previously, in [I], we considered the case when 

a'e(-oo,0)U(0,-^), (1.12) 

in which case the operator C has no discrete spectrum except zero. Under this condition we 
proved asymptotic stability for initial data close to a solitary wave both in the energy norm 
and in the weighted Banach norm, Li defined by, 



/ p 



[l + \x\ff{x)\\ LP . (1.13) 



In the present paper we extend this understanding to allow for the presence of additional 
discrete spectrum in the linearization: to be precise we will consider the case when 

, a a\/2(l + a/3) 

a e ( 7i^' 4c* >• (L14) 

In this case, there are, in addition to zero, 2 simple pure imaginary eigenvalues ±ifi, which 
satisfy the property 2/i > u>. If assumption fll.14}) is true for a fixed value uo, it also true for 

values of uj in a small interval centered at ujq. Let u(x,u) =1 J be the eigenvector of 

C associated to i/i. We choose the function u\(x) to be real, in which case U2{x) is purely 

imaginary. Then u* := I I , is the eigenvector associated to — i[i (see appendix A). We 

consider the initial value i^>(x, 0) = "00 (x) to be of the form 

ip Q (x) = e> e °^! UJo (x) +z u(x,uj ) + z u*(x,u ) + fo(x), (1.15) 



where /o belongs to the eigenspace associated to the continuous spectrum of C(luo). We assume 
that zq and /o are sufficiently small, and also assume a non-degeneracy condition which we 
now explain. Let (•, •) denote the Hermitian scalar product in L 2 of C 2 -valued functions, and 
(u, v) = u-\V\ + U2V2 for u,v G C 2 . Let E 2 [f, f] be the quadratic terms coming from the Taylor 
expansion of the nonlinearity: 

E*[f, f] = 5(x)[a'(C 2 )(f, f)*„ + 2a"(C 2 )(V^ /) 2 * w + 2a'(C 2 )(^^ /)/], / G C 2 . (1.16) 

The non-degeneracy condition has the form 

(E 2 [u{uo),u{uo)],T + {2ifio)) ^0, (1.17) 

where r+(2«/xo) is the eigenfunction associated to 2i/i = 2ifi(uo). This condition should be 
thought of as a nonlinear version of the Fermi Golden Rule of quantum mechanics ([HI Section 
XII. 6] or [5]), and will be referred to simply as the Fermi Golden Rule; see [TU| [TTj for the 
its introduction into nonlinear evolution equations. In appendix E we express ( 11.171) in terms 
of C and a(C 2 ), and hence show that the Fermi Golden Rule holds generically for polynomial 
nonlinearity. 

Our main theorem is following: 

Theorem 1.3. Let conditions [L~4\ ), ( TO)) and ( fOP hold, [3 > 2 andip{x,t) G C{R,H l ) be the 
solution to the equation U.fy) with initial value ipo(x) = ip(x,0) G H 1 fl Li of the form U.15\) 
which is close to a solitary wave e J 01 ^ W0 ." 

k(0)|<e 1/2 , ||/o||^<ce 3/2 . (1.18) 



Assume further that the spectral condition ftl.l4\ ) an d the non-degeneracy condition ([J.17| ) /io/d 



/or i/ie solitary wave with C = C(ujq) = Cq. Then for e sufficiently small the solution admits 
the following scattering asymptotics in C&(R) fl L 2 (]R): 

xj}{x,t) = e> v±(t) y w± {x) + e j ~ lLt $± + 0(t- u ), t -> ±00, (1.19) 

with some v > 0, where L = —-§^2, $± G C&(R) fl L 2 (M) are the corresponding asymptotic 
scattering states and ip±(t) = u±t + p-|-log(l + k±t) + Jt±, for some constants u±, p±, k±, x±. 

The asymptotics (11.191) can be rewritten in terms of the original complex notation as: 

i){x,t) = e lLp±{t) i) u)± {x) + W{t)(l) ± + 0{r v ), *->±oo, (1.20) 

where W(t) is the dynamical group of the free Schrodinger equation, and <f>± = $^_ + i<&\ ($±, 
k = 1, 2, being the components of the vector-function $±). Thus the main conclusion is that 
asymptotically the dynamics decomposes into a nonlinear bound state e tip± ^ip UJ± (undergoing 
uniform phase rotation, modulo the logarithmic phase shift in (p±(t)) and a solution of the free 
Schrodiniger equation W{t)(j>±. 

The proof is divided into the following three main steps which are carried out, respectively, 
inii and §5J 

Step 1 To decompose the solution ij;(t) according to the spectral decomposition of the operator 
C given in §2.1[ and obtain a system of equations in §3. II equivalent to (II. ip . This system 
is then put into a canonical form, in which certain non-resonant terms are excluded. This 
is carried out in §3.41 the final form of the equations being given in §3.4.51 



Step 2 To use the time decay for the linearized evolution on the continuous spectral subspace 
established in §2.21 to prove asymptotic stability of the solitary waves in §4.61 



Step 3 Finally, to construct the wave operator and obtain the scattering asymptotics (11.191) . in 



2 Linearization 

In this section we summarize the spectral properties of the operator C which appears in the 
linearization of ( 11. 21) about a soliton, and then give some estimates for the linearized evolution. 
The proof of these properties can be found in pQ, with the exception of proposition 12.31 which 
is proved in appendix C. 

2.1 Spectral properties of the linearization 

The linearized equation reads 



X(x,t) = Cx{x,t), 



c==r'B = (_ D « D l). (2.1; 



Theorem 11.11 generalizes to the equation (12. ip : the equation admits unique solution x( x ^) £ 
Cft(R, H 1 ) for every initial function %(x, 0) = Xo £ H 1 . Denote by e ct the dynamical group of 
equation (12. ip acting in the space H 1 ; for T > there exists Ct > such that 

||e ct Xo||^<c r ||xo||Hi, \t\<T. (2.2) 

The resolvent R(A) := (C — A)" 1 is an integral operator with matrix valued integral kernel 

R(A, x, y) = T(A, x, y) + P(X, x, y), (2.3) 



where 



/ \ / e ik+\x-y\ _ e ik+(\x\ + \y\) _x( e ik+\x-y\ _ e ik+(\x\+\y\)\ 

T(X,x,y) 

e ik-\x-y\ _ e ifc_(|x|+|y|) j/gifc_|x-y| _ e ik-(\x\+\y\)\ 



\ 4fcl AkI J 



(2.4) 



1 / e ik +W e ik -W \fia-2k- i(3 \ f e ik +^ -ie ik +^ 

n\x,y) ~2Dy ie ik+\x\ _ ie ik.\x\ J y _ %f3 -ia + 2k + J \ e ik ~ M ie ik ~ M 

• (2.5) 

As explained already in §1.21 the continuous spectrum consists of C + U C_ , and correspondingly 
&±(A) = V~ w T iX is (respectively) the square root defined with a cut in the complex A plane 
so that fc±(A) is analytic on C\C± and Imfc±(A) > for A G C \ C±. The constants a, f3 and 
D = -D(A) are given by the formulas 

a = a + a'C 2 , (3 = a'C 2 , D = 2ia(k+ + fc_) - Ak+k. + a 2 - (3 2 . 

In addition to this continuous spectrum, there is discrete spectrum, which appears in this 
formalism as the set of poles of R(A); these poles in turn correspond to the roots of the 
determinant D(X). From the analysis in [T] we know that if a' G (a/v^C 2 , a/C 2 ), then the 



determinant has the following roots: A = with the multiplicity 2 and two pure imaginary 
roots 

A = ±i/j, = ±i Wa 2 - P 2 , fi<u. (2.6) 

Note that spectral condition (j!.14p is more restrictive. It implies in addition that 2/i > u, as 
can be verified by a simple computation. 

The generalised null space X° of the non-self-adjoint operator C is two dimensional, is 
spanned by j^, d^^, and 

CjV u = Cdjl „=&„,. 

The symplectic form Q for the vectors ip and rj is defined by 

Sl(il>,Ti) = (il>,jri) (2.7) 

By Lemma [1.21 

n(j*u>, d^J) = -\d u J \^\ 2 dx ± 0. (2.8) 

Hence, the symplectic form Q is nondegenerate on X°, i.e. X° is a symplectic subspace. There 
exists a symplectic projection operator P° from L 2 (M) onto X° represented by the formula 

P °^ = /* l* s\ U> >&»*<*)&» + W> *J)d u * u ]. (2.9) 

Remark 2.1. Since j^, cL^w ^e in /7 1 (IR) i/ie operator P° extends uniquely to define a 
continuous linear map H~ l (M.) — > X°. In particular this operator can be applied to the Dirac 
measure S(x). 

Denote by X 1 the eigensubspace corresponding to eigenvectors u and u*, and by P 1 a 
symplectic projection operator from L 2 (R) onto X 1 . It may be represented by the formula 

rV = ^4„+«^L. (2 . 10) 

{u,ju) (u*,ju*) 



since (u,ju*) = 0, and (u,ju) = (u*,ju*) ^ by (I3.66p . Finally, P c = 1 — P° — P 1 is the 
symplectic projector onto the continuous spectral subspace. 

2.2 Estimates for the linearized evolution 

We now recall from [1] some estimates on the group e ct which will be needed in §H First we 
recall a bound for the action of e ct on the Dirac distribution 5 = 5(x) for small t. Lemma 8.1 
from [I] gives the following small t behaviour: 

||e ct 5|Uco = £>(r 1/2 ), t — 0. (2.11) 

Second we list the large time dispersive estimates. To do this let us introduce, for (3 > 2, a 
Banach space M.p, which is the subset of distributions which are linear combinations of L\ 
functions and multiples of the Dirac distribution at the origin with the norm: 

||^ + CT(x)||^:=||^|| L i + |C|, (2.12) 

and a Banach space Bp = B(Ai f3 , L ?*) as the space of continuous linear maps M.p — > L% for 
any (3 > 2. 



Proposition 2.2. (see ITj) Assume that the spectral condition ftl.l4\ ) holds. Then for h E Mp 
with j3 > 2 the following bounds hold: 

\\e ct P c h\\ L ^<c(l + t)- 3 / 2 \\h\\ M0 , (2.13) 

lle^ir^lU- < c(l +t)- 3 ' 2 \\h\\ Mf> , (2.14) 

\\ e ct c -i n ± h ^ < c (i + t )-W\\h\\ Mp , (2 15) 

where II + (resp. TI~) is the spectral projection operator onto the spectral subspace associated 
to C + (resp. C-), the positive (resp. negative) part of the continuous spectrum. 



We shall also need the following bound, which is new. 



Proposition 2.3. Assume that the spectral condition jjjl.l4\ ) holds. Then for h E M.p with 
(3 > 2 and t > 1 the following bounds hold: 

\\e ct (C T 2ifM - O)- l P c h\\ L ~ < c(l + t)-^ 2 \\h\\ Mp , (2.16) 

||e C4 (C + 2ifi - 0)- 1 H ± /i.|| x ,oo g < c(l + t)- 3/2 \\h\\ Mp . (2.17) 

We prove this proposition in appendix C. 

3 Spectral decomposition and canonical forms 

In this section we will use the spectral resolution associated to the operator C to decompose the 
solution ip, obtaining evolution equations for the different spectral components. Then, following 
[2], we introduce normal coordinate transformations to transform the evolution equations into 
simpler canonical forms in which certain non-resonant terms are absent. The final form of these 
equations is given in §3.4.51 

3.1 Modulation equations 

In this section we present the modulation equations which allow a construction of solutions 
ip(x, t) of the equation ( 11.11) close at each time t to a soliton i.e. to one of the functions 

in the set S described in §1.21 with time varying ("modulating") parameters (cu, 6) = (u(t), 9(t)). 
We look for a solution to (jl.2p in the form 

il>(x,t) = e? e ®$(x,t), $(x,t) = ^ u{t) (x)+x(x,t). (3.1) 

Since this is a solution of (11.21) as long as x — with 9 = uo and to = it is natural to look for 
solutions in which x is small and 

0(f) = / cu(s)ds + j(t) 
Jo 

with 7 treated perturbatively. We look for x = w (a ; , t) + f(x, t) where w = zu + zu* E X 1 and 
/ E X c . Now we give a system of coupled equations for u(t), j(t), z(t) and f(x,t). 



Lemma 3.1. Given a solution of (11.21) in the form (13. ip with f G X c as just described, the 
functions uj{t), ■yit), z{t) and f(x,t) satisfy the system 

a (P°Q[x],g) f32) 

<9 W * W - flLP°x, *> ' l j 

. = (Q[xU(aA-3,P°x)) r „ - 

7 {d u * u - flLP°x, *) ' l ' j 

(uju)(z-i/j,z) = (Q[x],ju) - (d^ - d^P 1 f , ju)u - (x,u)7 ( 3 - 4 ) 

/ = C/ + c^P c x - 7P C (JX) + P C QM, (3-5) 

where Q[%] = — 5(x)j~ 1 (F(\l/ a; + x) — F(\E f a; ) — F'( 1 S r £J )x) represents the nonlinear part of the 
interaction. 

Proof. This can be proved as Proposition 2.2 in [2]. □ 

3.2 Frozen spectral decomposition 

The linear part of the evolution equation (J3.5P for / is non-autonomous, due to the dependence 
of the operator C on u>(t). In order to make use of the dispersive properties obtained in §2.21 
it is convenient (following [2]) to introduce a small modification of (13. 5p . which leads to an 
autonomous linearized equation. Let us fix an interval [0, T] and decompose /(£) G X£ into the 
sum 

f = g + h, ff e4 hex?.. (3.6) 

Here X^ = P^X is the spectral space associated to the discrete spectrum at time T and Xf. = 
P^X is the spectral space associated to the continuous spectrum at time T, P^ = P c (u(T)) 
and P^ = I — Py. In the following, we denote ujt = w(T) and Cr = C(u>t)- We will obtain 
estimates uniform in T, and later consider the limit T — ► +oo. 

We now introduce a shorthand for the bounds we are about to prove: 1Z(A, B, . . . ) (resp. 
1Z(u, A, . . . )) is a general notation for a positive function which remain bounded as A, B, . . . 
approach zero (resp. if u is close to u and A, . . . approach zero); it could be unbounded and 
even infinite if u is outside some vicinity of ujq. The formula / = IZg implies that |/| < IZg. 
Introducing also the notation 1Z\(uj) = 7Z(\\uj — uj \\c[o,t}), we get 

Lemma 3.2. The function g is estimated in terms of h as follows: 

\\g\\L™ = TIi(u>)\uj - u>T\\\h\\ L ™ (3.7) 

Proof. Let us use the identities 

P d (g + h) = 0, P d T g = g, P d T h = 0. 

Then we get 

g + (P d - P d T )g + {P d - P d T )h = 0, 

and P d — Pf. is a "small" finite dimensional operator: 

\P d - P d T \ = \P d (u t ) - P d (u T )\ < max w . 6(W)Wr) |9 w P d (a;*)||a; - u T \. 



a 



Applying the projection P^ to (I3.5I) . we get 

h = C T h + P£,[(C - C T )/ + P C Q[ X ] + cod u P c X - iP c (jx)}- (3.S 

8 



3.3 Asymptotic expansion of dynamics 



The preceding sections have provided a change of variables ip i— > (u),^,z,h) under which (11.21) 
is mapped into the system comprising (I3.2l) . (l3.3p . (13.41) and (13.81) . Since we are interested in 
proving that for large times z, h are small it is necessary to expand the inhomogeneous terms in 
these equations in terms of z, h. This is carried out in this section, leading to the conclusion that 
the system (I3.2p . (l3.3p . (13. 4p and (13.81) can be written in more detail as the system comprising 
(E22]),(J329D, (HLH and <^M>- 



3.3.1 Preliminaries 

This section is devoted to some useful preliminary estimates. We start with a bound for the 
denominator (d u ^f — d LU P°x, $), where \l/ = ty^, that appears in the equation of motion (13. 2h - 
( 1331) . We have, with A = (d^, *), 



(d u V-d u P x,$) = (d u V,V)(l 

with 



(d^, X )-(d UJ P x,<!>) 



(d„*,V) 



A 1 



(d^, X )-(d UJ P x^) 



(d^, X )-(d^P°x,$) 



K{u) (\z\ + H/IU- + ||/||£„ 



(3.9) 



(3.10) 



A s ' v ' "" "~-p "" ""-?, 

We also need to expand the nonlinear term F(^) = a(\i(j\ 2 )ip near the solitary wave since the 
inhomogeneous terms all involve E[x], the nonlinear part of 5(x)F, defined using the Taylor 
expansion of 5(x)Fif) near \1/: 



(J(x)F(V) = 5(x) a(C 2 )V + a(C 2 ) X + 2a'(C 2 ) ( X , ¥)* + E[ X ]. 



(3.H) 



Thus E[x] contains all the higher order terms which are at least quadratic in x-, as X ~~ * 0, and 
Q[x] = J-^M- We expand -E[x] i n the form 



E[x] = E 2 + E 3 + E_ 



R, 



where Ej is of order j in x an d E^ the remainder. It is easy to check that 

E 2 [x,x}=S(x)\a'(C 2 )\x\ 2 ^ + 2a"(C 2 )^,x) 2 ^ + 2a\C 2 )^, X )x 



(3.12) 



(3.13) 



^ 3 [x,X,x]=5(x)[a'(C' 2 )|x|\ + 2a // (C (2 )(vl/, x )\ + 2a / '(C (2 )(v[/, x )| x | 2 v[/ + -a w (C' 2 )(vl;, x ) 3 ^ 

(3.14) 
For En we have 

^ = ^(c,|^|,|/(0)|)(|^| 4 +|/(0)| 4 ), (3.15) 

It also useful to define E 2 [xi, X2], (resp. E 3 [xi, X2, Xs\) as a symmetric bilinear (resp. trilinear) 
form 



E 2 [xi, X2] = 5(x) [a'(C 2 )( Xl , X2)* + 2a" (C 2 )^, X i)(% X2)* + a'(C 2 ) ({% X2 ) X i + (% Xi)x 2 

(3.16 

E 3 [ X i,X2,Xz] = *(a;)[io'(|*| 2 )X;(xi,^)Xfc + ^o"(|*| 2 )X;(W, Xi )(* s Xi)Xfc (3-17) 

+ V(|vI/| 2 )^(xl/, x ,)( Xj , Xfc )vl/ + V(|^| 2 )(^Xi)(^X2)(^X3)^ 



Here summation is taken over all transposition of integers 1, 2, 3. Notice also that 

(E 2 [X,Y],Z) = (X,E 2 [Y*,Z}) (3.18) 

where X, Y, Z, are complex valued vector functions and Y* = (Y\, Y 2 ). 

In the remaining part of the paper we shall prove the following asymptotics: 



-0 



r 1 , z(t)~r 1/2 , ||w(t)||^i ~r 1/2 , t^oo. (3.19) 



Remark 3.3. To justufy these asymptotics, we will separate leading terms and remainders in 
right hand side of equations l[3.ty) - (3~4\ ), Ii3.8\) . Namely, we shall expand the expressions for Co, 
7 and z up to and including terms of the order 0(t~ 3 ^ 2 ), and for h up to 0{t~ x ) keeping in mind 
the asymptotics h'J.19\) . This choice is necessary for applicaton of the method of majorants. 

3.3.2 The equation for uo 

Using the equality Q[x] = jE[x], and the fact that j(P )* = P°j (where * means adjoint with 
respect to the Hermitian inner product ( ■ , ■ )), we rewrite 

(P°Q[X],$) = (P°jE\x],$) = -(E\x],j{P°)*$) = -<£[*], PW 

with x — w + / an d $ = v I/ + x = ^ / + w + /- Then equation (13.21) for Co can be expanded up 
to C(t~ 3 / 2 ), assuming (13.191) . as follows: 



1 



(E 2 [w,w] + 2E 2 [w,f] +E 3 [w,w,w],jW) + (E 2 [w,w],P°jw) 



(3.20) 



1 



2 (£ 2 [w,w],j^)((^,w>-(^P u w,^>)+^ 



where 

n R = n(u, \z\ + ||/|U~ )(M 2 + ||/|U~ ) 2 (3.2i) 

Substituting w = zu + zu*, we can write (I3.20p in the form 

Co = tt 20 z 2 + tt u zz + tt 02 z 2 + tt 30 z 3 + tt 21 z 2 z + tt l2 zz 2 + tt 03 z 3 + z(f, tt' w ) + z(f, Q,' 01 )+n R (3.22) 

Let us now display explicitly some important terms of this expansion. First we compute the 
quadratic terms in (13.201) which are of order t _1 according to (13.191) . They are obtained from 
the term 

(E 2 [w,w],jV) = z 2 {E 2 [u,u),j^} + z 2 (E 2 [u*,u*]jy) + 2zz(E 2 [u*,u]jy). (3.23) 

in expression ( 13.21) . Let us take into account definition of E 2 , the identity (\I/, j\I/) = and that 
$ = (0,0), w = zu + ~zu* with u = (ui,u 2 ) where u\ real and u 2 pure imaginary. Then we 
obtain 

(E 2 [w,w]J*) = (5(x)2a'(C 2 )^,w)w,j^) = 2a / (C 2 )(z + z)(u(0)^(0))(z-^(u(0), 3 ^(0)) 
= 2(z 2 -z 2 )a\C 2 )(u(0)MQ))HO),mO)). 



Therefore 



n 20 = Qo2 = JE2[u,u],j*) = _|. o / (c .2 )(u(0) ^ (0))(u(0)) ^ (0)) (324) 
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is purely imaginary and 



nu = _ 2 <g.i".ni*> =0 . 



A 



Using the property (13.181) . we find that 



E 2 [u*,jV] 



0' — o' — —9 

is 10 — ii 01 — z 



(3.25) 



(3.26) 



Remark 3.4. Since f e Xf £/jen 



(/,nio> = <P c /,n / io> = </,j'P c r 1 n / io> 



Therefore we can substitute Q' 10 in V3.2ty by their projection jP c j 1 / 



io- 



3.3.3 The equation for 7 

Using again the equality Q = j ] E we get 

(Q\x]J{d u *-d u P X )) = (E\x],d u *-d u P°x) 
Therefore (ET3J1 . (13T91 . (I3~T0|) . (15321) . ([335]) imply 
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A- 1 



(E 2 [w,w) + 2E 2 [w, f] + E 3 [w,w,w},d^) - (E 2 [w,w},d UJ P w) 



(3.27) 



- A- 2 ( J E 2 [w,w],^^)((^^,w)-(^P°w,^))+r i? , 
where 

T R = K(u, |*| + |/(0)|)(M 2 + |/(0)|) 2 = fc(u,, |*| + ||/|U-,)(k| a + ll/IU^) 2 (3.28) 

since |/(0)| < ||/||l~ • Equation (13.271) can thus be represented in the form 

7 = r 2 o* 2 + r 11 ^+r # 2 + r 3 o2 3 + r 21 ^+r 12 ^ 2 + ro32 3 + ^(/,r / 10 ) + 2(/,r , 01 ) + r i? , (3.29) 

where 



(E 2 [u,u],d„V) r _ n (E 2 [u,u%d^) (E 2 [u*,u%d w V) 

1 20 — a j *■ n — z , 1 02 



A ' xx A 

1 in — -^ . 1 -L m — ^ 



A 



(3.30) 
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') A 01 



A ' U1 A 

3.3.4 The equation for z 

Denote x = (u,ju) and rewrite (13.41) in the form: 

(E 2 [w,w\ + 2E 2 [w, f] + E 3 [w,w,w],u) 



Z — I/jLZ 



X 



(3.31; 



(d u wJu){Eh[w,w],jV)-{w,u){E2[w,w],d u V) 

+ xA + Zr > 
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where 



Z R = K(u;,\z\ + \\f\\ L ~)(\z\ 2 +\\f\\ L ~y 



Equation f 1 3 . 3 1 j) can be represented in the form 

z = ifiz + Z 20 z 2 + Z n zz + Z 02 z 2 + Z 30 z 3 + Z 21 z 2 z + Z 12 zz 2 + Zq-^zT 
+ z(f,Z' 10 )+z(f,Z' 01 ) + Z R , 

where, using the calculations in the previous two sections, we have in particular, 



Z 2 o 
Z21 



(E 2 [u,u],u) 



X 



'11 



\E 2 [u,u%u) 



x 



^ 02 



(E 2 [u*,u*],u) 



x 



(3E 3 [u*,u,u],u) 



x 



(3.32) 



(3.33) 



+ 



(duU*,ju)(E 2 [u,ulJ^) ~ (u*,u)(E 2 [u,u},d^) - (u,u)(2E 2 [u*,u},d u> ^) 



xA 



Z[ 
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2 E 2 [u\u] = E 2 [u,u] 



x ' ~ 01 



X 



(3.34) 



3.3.5 The equation for h 

We now turn to equation ( 13.81) for h that we rewrite in the form 

h = C T h + P^ [(C - Cr)/ + P c jE 2 [w, w] + 7P C J-V + H R 

where remainder H R is 

H R = ¥ c j{E[ X ]- E 2 [^M) +ud ul Y c x + i 1 P c r l ™ 

= P c j(E[ X ] - E 2 [w, w]) - ud^x + 7i _1 w - 7 P d J^w 

For the .£/# we have, recalling (12.121) . the following estimate 



\H 



RWMp -y.~-.V-X ■ II- l.-.fl/M", ■ ,"...- ..-_„ „. „-_ . 



lZ(cu)\j\\z\ = lZ(u>, \z\ 



L~ g K\zr + \z\u\\ LT0 + 



(3.35) 



(3.36) 



K(u>, \z\ + ll/IU- )(|z| 3 + kill/Ik-, + H/lli- ) + W(«)|a;|(W + ||/||j/-,) 



(3.37) 



Now we continue the isolation the leading terms in the right hand site of (I3.35p . Note that, 
from the formulae in the discussion surrounding (II. lip . 



C-C T = r\uj-ujT)+r 1 (V -V T ), where V = -5{x)[a + bP x ] 



Also 



P^P C = Pr[Pr + Pr - P d l 



c i-rtd 



Pr + Pr[P 



Tl 1 - T 



Therefore f]3.35f) becomes 



h = C T h + a{t)¥ c T r l h + P c T jE 2 [w, w] + Hj 



R 



(3.38) 



with a{t) = oj — uj t + 7 , and 

H' R = P C T [H R + a(t) 3 - l g + j~ l (V - V T )f + (P d T - P d )j(E 2 [w, w] + 7/)]. 
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Using the identity P^ = 1 — P!f, we obtain 
WrWm, < fci(w, 1*1 + ll/IU^)(k| 3 + klll/IU-, + H/lli- + |w - ta*|(|*| a + |7|||/IU^). (3.39) 

Next we need an additional construction to combine first two terms in RHS of (13.381) . Namely, 
lemma [3751 below shows that the "main part" of the second term is ia(t) (Ilj — II^)/i, where II + 
and n~ are the spectral projection operators on the spectral space associated to the positive and 
negative part of the continuous spectrum respectively at time T; see the discussion preceding 
(j2.14p . Hence, we denote 

C M (t) = C T + ia(t)(Il+ - n~) (3.40) 



and rewrite (13.381) as 
where 



h = C M (t)h + P c T jE 2 [w, w] + H R (3.41) 



H R = H' R + o(t)\V c T r l - ^(n+ - U T )]h (3.42) 



Lemma 3.5. For h £ X% we have 

||[p^- 1 -z(n+-n T )]/ i || L i < \\h\\ L ^. (3.43) 

This lemma is proved in appendix D. Lemma [3.51 and the bound (13.391) imply 
Proposition 3.6. The remainder H R admits the bound 
\\H R \\ Mp < Mu, \z\ + ll/IU- )(|z| 3 + \z\\\f\\ L oo + ll/Hi. + \u - u T \(\z\ 2 + ll/IU- )). (3.44) 

3.4 Canonical form of the equations 

Our goal is to transform the evolution equations for (u, 7, z, h) to a more simple, canonical form. 
We will use the idea of normal coordinates, trying to keep unchanged the estimates for the 
remainders as much as is possible. This means we observe that for our purposes the unknowns 
(lj, 7, z, h) lie in a neighbourhood of the point (o>o, 0, 0, 0) in the space GlxlxCx Z/^g. We 
seek a change of variables 

: (w, 7, z, h) ^ (wx, 71, Zi, h{) 

such that is a diffeomorphism between neighbourhoods of (ujq, 0, 0, 0) in the space G M x 
M. x C x L%, and Z)6(wo,0,0,0) is the identity. This map is obtained, as usual in the 
normal form method, by looking for (o>i, 71, z%, hi) as a power series in (u, 7, z, h), starting with 
the identity map at highest order - see (137471) - (EJ749J) . (13371) . flSZZSD and (137821) for the detailed 
expressions. The coefficients in these expressions are then chosen so as to put the equations 
for (u>i,7i, zi, hi) in a simpler canonical form which is suitable for our subsequent work. This 
final canonical form for the system is summarised in section 13.4.51 

3.4.1 Canonical form of the equation for h 

As a starting point we expand out the middle term on the right hand side of (I3.4ip . obtaining 

h = C M (t)h + H 20 z 2 + H n zz + H m z 2 + H R . (3.45) 
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Here, the coefficients Hij are denned by 

H 20 = P c T jE 2 [u,u], H n = 2F c T jE 2 [u,u*], H l2 = P c T jE 2 [u*,u*]. (3.46) 

We want to extract from h the term of order z 2 ~ t~ l . For this purpose we expand h as 

h = h 1 + k + k 1 , (3.47) 

where 

k = a 20 z 2 + auzz + a 02 ~z 2 , (3.48) 

with some a^ = %• (o>, x) satisfying a^ = <%, and 



fci = -expf fc M (r)dr)k(0). 



(3.49) 



Note that ki is just the solution of the corresponding homogeneous equation ki = Cj^ki, since 
the operators Cm(^) all commute for different values of t. It follows from &i(0) = — k(0) that 
^(0) = h(0). 



Lemma 3.7. There exist aij e L^a in l[3.48\ ) such that the equation for h\ has the ft 



form 



h = C M {t)h! + H R (3.50) 

where Hr = Hr + H' , with estimates as in (j3.44p . and also 



» =fti(w)|*| 2 . (3.51) 



Proof, (cf. Section 4.2.2 in [2]) We substitute (I3.48P into (j3.4ip and equate the coefficients of 
the quadratic powers of z. In addition we replace the discrete eigenvalue /x(t) by its value at 
time T, i.e. \xt = //(o;(T)), and include the correction in the remainder. Then we get 



20 ~ 2in T a 20 - 


- — CtO20 


H u = 


= —Cyan 


02 + 2ifi T a 02 = 


= — CtO02 



and 

H R = H R + H 

where H' is defined as 

H' = ^duOvKiu, \z\ + \\f\\ L ^)\z\ 2 (\z\ + \\f\\ L ~ ) 2 (3.52) 

+ ^di/Hfa \z\ + ||/|| L » )\z\{\z\ + ll/IU- ) 2 + ^a^Ml^lVr - /<| - ^(11+ - H^)fc. 

The dependency in x appears here through the coefficients a^ = a,ij(uj,x). Notice, from (11.161) 
and the formuae for the projection operators in §2.11 that each H^ e X? is the sum of a multiple 
of 5(x) and a function exponentially decreasing at infinity. Hence, there exists a solution an in 
the form 

a u = -C^Hn, (3.53) 
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where C^ 1 stands for regular part of the resolvent R(X) at A = since the singular part of 
R(X)Hn vanishes for H\\ G X£. The function an is exponentially decreasing at infinity. 
For a2o and Oq 2 we choose the following inverse operators: 

a 20 = -(Or - 2tfi T - 0)~ l H 20 , a 02 = a 20 = -(C T + 2i// T - O) -1 ^, (3.54) 

This choice is motivated by proposition 12. 3[ and putting t = in that proposition we have the 
bound (13.511) . The remainder H 1 can be written as 

H' = ^(C T - 2i/ptm - Q)~ l A m , me {-1,0,1} (3.55) 

with A m G Xf., satisfying the estimate 

\\A m \\ Ll = K(u, \z\ + ||/|| L ~ )|*| (\z\\u T -u\ + (\z\ + ||/|U~ ) 2 ). (3.56) 

D 

3.4.2 Canonical form of the equation for u> 

We want to remove all terms in the right hand side of (13.221) except the remainder Q R . This is 
possible by methods of Buslaev and Sulem [21 Proposition 4.1] since Qu = by ( 13.251) . 

Lemma 3.8. There exist coefficients bij(ou), < i,j < 3, and vector function b^Ax^u), < 
i,j < 1; such that real-valued function uj\ defined as 

ui = lo + b 20 z 2 + b 02 ~z 2 + b 30 z 3 + b 21 z 2 l + b 12 zz 2 + b m ~f + z(f, b' w ) + z(f, b f 01 ) (3.57) 

obeys a differential equation of the form 

vi = Ci R , (3.58) 

where Qr satisfies the same estimate $3.21\) as Qr: 

Q R = K(ou, \z\ + ||/|U- )(|*| a + II/IIl^) 2 (3.59) 

Proof. The calculation follows the classical method of normal coordinates. Substituting u, z, 
and / from (13.221) . (I3.33p . (I3.5P into the equation for Cj\ and comparing the coefficients of z 2 , 
zf, etc. leads to a system of equations for the coefficients 620, b' w , ets. (cf. [21 Proposition 4.1]) 

tt 20 + 2^620 = 0, (3.60) 

n' w + iiib' w + C*b' w = 0, 

n 21 + 2Z 11 b 20 + ifj,b 21 + 2Z 20 b 02 + (F n , 6i ) + (F 20 , b' 01 ) = 0, 

tt 30 + 2Z 20 b 20 + Sifjbgo + (-F20, fr'io) = °- 

It should be noted that the resonant term zz in the equation for Co\ is absent. From the first 
equation of (I3.60p we obtain 

b 20 = b 02 = —Q 20 , (3.61) 
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Multiply the second equation of (I3.60P by j we get 

jQ' 10 + i/jjb' 10 - Cjb' 10 = 

since jC* = —Cj. Without loss of generality we can assume that jT2' 10 G X£ by Remark 13.41 
Therefore, there exists the solution b' w in the form 

&io = Vox = ~j(C - t^y 1 ]^ (3.62) 

where (C — i/i) -1 stands for regular part of the resolvent R(\) at A = ifi since the singular part 
of R(X)jQ' 1Q vanishes for jT2' 10 G X^. The functions 6' 10 , b' m decrease exponentially at infinity, 
and the equations for b 2 \ = bi 2 , 630 = &03 can be easily solved. □ 

3.4.3 Canonical form of the equation for z 

In this section we obtain a canonical form of the equation ( 13.331) for z, and carry out a com- 
putation of the coefficient of the resonant u z 2 ~z " term, which gives the Fermi Golden Rule. 
Substituting (13.61) and (13.471) into (13.331) and putting the contribution of g + hi + k\ in the 
remainder Zr, we obtain 



z = ifiz + Z 2 oz 2 + Z u zz + Z 02 z 2 + Z 30 z 3 + Z 21 z 2 z + Zi 2 zz 2 + Z 03 z 3 (3.63) 

+ Z^qZ + Z 2 \Z z + Zy^zz + Zq^z + Zr. 

We have by (13T47D - ( 13^481) 

Z' 30 = (a 20 ,Z' w ), Z' 21 = (a n ,Z' 10 ) + (a 20 ,Z' 01 ), (3.64) 

-^03 = ( a 02,Z 0l ), Z 12 = (a 02 ,Z 10 ) + (a n ,Z 0l ). 

We are particularly interested in the resonant term Z' 21 z 2 ~z. Formulas ( 13.341) . ( 13.461) . ( 13. 53ft . 
(ET53J) imply 

Z' = -(Cz 1 2P c T jE 2 [u,u% 2 E2 ^ U *h - ((C r -2^ T -0)~ 1 P^E 2 b,d,2 E2 ^^ ) (3.65) 

X X 

For the coefficient x = k[lS) we get (see 0, Proposition 3.1]) 

x = (u,ju) = i5, with <5 > 0. (3.66) 

Now we can prove 

Lemma 3.9. Suppose that the non- degeneracy condition \1.11\) is satisfied, then 

ReZ 21 <0 (3.67) 

for uj in some vicinity of ujq. 

Proof. We first notice that the coefficient (C^ 1 2P^jE 2 [u,u*), E 2 [u,u*}) appearing in the ex- 
pression ( I3.65P for Z' 2l is real, since the operator C^ 1 2P^j is selfadjoint. Hence by ( 13.661) 
ReZ 21 reduces to 



ReZ , Re2 ((C r ~ 2tfi T - 0)- 1 P c T jE 2 [u, u], E 2 [u, u]) 



x 

--Im {R(2in + 0)P c T jE 2 [u, u], E 2 [u, u]), 
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where we denote 

R(\)=R T (\) = (C T - \y\ Re\>0, and [i = [i T - 

Using that P§, commutes with R(2ifj l + 0), we have R(2i(j, + 0)P% = P c T R{2ifi + 0)P c T . We have 
also that (Py)* = —jP^j, hence 

2 
ReZ' 21 = -lm(R(2i(j, + 0)a,ja) 

with a = Pj.jE2[u,u]. The function A h- > (R(X)a,ja) is analytic in the region C \ (C+ U C_) 
since a G X£. Hence by the Cauchy residue theorem we have 

W2, A , + 0)a,ja) = -— J dX ^ /.^V l (3-68) 

c+uC- 

Now we use the representation 

«(A + 0)-K(A-0) = - ^'y - fc(> ^ W J, Ae C± , (3.69) 

where D = L>(A + 0) and fc+ = A;+(A + 0) < for A G C+, fc_ = A;_(A + 0) > for A G C_ 

r± (A) = (De tk±lxl - De~ ik± W)v ± + 4(3ik±e ik ^v T , s ± {X) = sm(k±x)v± (3.70) 

are the even and the odd eigenfunctions of the operator Cy corresponding to the point A G C± 
of the continuous spectrum (see appendix B). The representation (13.691) can be checked by 
direct calculation using formulas (ID. 2[) - flD. 4D for the resolvent. Then equation (13.681) becomes 
(with A = iv) 



— UJ 

<J«2fc + 0)<»,,«> = lar J k _ ]D? 


(r- 


,ja){r-,ja) 

v — 2/i 


oo 

167T y 


— OO 






U) 


since the function a is even. Using that 








1 




1 „, - 





dv (T + ,ja)(r + ,ja} 
k + \D\ 2 ^-2/i + iO 



v + iO v 

where p. v. is the Cauchy principal value, we have 



p.v. iix5{v) 



oo 



(R(2ifi + 0)a,ja) = -—- / - — — T7^P- V - 



16vr J k_\D\ 2 v-2n IGn'J k + \D\ 2 v - 2/i 



i (T+(2ifj,),ja)(T + (2ifi),ja) 
16 k + (2ifi + 0)\D{2ifi + 0)\ 2 l ' ' 

The integral terms in (13.711) is real. Thus, 

|(r + (2^ r ),£ 2 [u,u]>| 2 



Im (R T (2ifi T + 0)a,ja) 



lQk + (2i{i T + 0)\D(2ifi T 



The non-degeneracy condition (11.171) implies that (t + (2z'//t), E 2 [u, u]) ^ in some vicinity of 
ujq. Using also the inequality k + (2ifi T + 0) < 0, we deduce ReZ 21 < 0. D 
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We now need an estimate on the remainder Zr. 
Lemma 3.10. The remainder Zr has the form 

z R = n 1 (cu,\z ] ■ 



l? ) (M +11/11^) +I 2 II u; t'-^|||/i|U^ 3 + |2;|||A;i||loo /3 + |z|||/ii||l^ 3 



(3.72) 



Proof. The remainder Zr is given by 

Z R = Z R + z(f - k, Z' w ) + z(f - k, Z' m ) 
where Z R satisfies estimate (13.321) . Since / — k — g + k\ + hi, we have by (13.71) 



\(f-k,Z[ )\ < H(u>){\\g\\v-+ Wkth- +||/ii|U-,) 



< TZi(u)(\ut — wlll/il 



l*i| 



Mb 



'-0 ' " " ,|iJ -/3 

which implies (13.721 

We can apply now the method of normal coordinates to equation (13.631) . 

Lemma 3.11. (cf. ]H Proposition 4-9]) 

There exist coefficients Cij such that the new function z\ defined by 

z\ = z + c 20 z 2 + cuzz + c 02 ~z 2 + c 30 z 3 + Ci 2 ZZ 2 + c 03 z 3 , 

satisfies an equation of the form 

zi = ifj,(u)zi + iK(u)\zi 1 2 zi + Z R 

where Z R satisfies estimates of the same type as Zr, and 

Re iK = ReZ' 21 < 0. 



□ 



(3.73) 



(3.74) 



(3.75) 



Proof. Substituting z\ in equation (13.631) for z and equating the coefficients, we get, in partic- 
ular, 



and 



1 7 % 7 l 7 

C20 — — ^20) c ll — ^11; c 02 — — 7^ z '02 

\X /i 3/X 



iK — Z 2 \ + Z 2 i + C\\Z 2 q + 2c 2 qZh + C\\Z\\ + 2Cq2^ 



02 



(3.76) 



(3.77) 



It is easy to check that all of the coefficients Zn, Z 2 o, Z 02 , and Z 2 ± defined in (13.341) are pure 
imaginary, and hence (I3.75P follows immediately. □ 

Denoting Kt = K(u>t), the equation for z\ is rewritten as 



where 



Zr + TZ^u, \z\ 



Z\ = ijjLZi + iK T \zi\ 2 zi + Z 



l™)\z\ \ut -uj\= TZi(v, \z\ 



r (3.78) 



+ ii/iu-; (n 2 + ii/iu-; 2 (3.79) 



+ \z\\ojt — o;|(||/i||loc , g + \z\ ) + |z|||/ci||l^° + |^|||/ii||z,° 
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It is easier to deal with y = \zi\ 2 , rather than z\, because y decreases at infinity while z\ 
is oscillating The equation satisfied by y is simply obtained by multiplying (I3.74p by z\ and 
taking the real part: 

y = 2Re (iK T )y 2 + Y R , (3.80) 

where 

Y R = n^uj, \z\ + \\f\\ L ~)\z\ \(\z\ 2 +\\f\\ L ~f+\z\\uj T -uj\(\\h\\ L ^ + \z\ 2 ) + \z\\\k 1 \\ LT ■ +|*||MUo 



-P' 



-0' "-"-J3 '-""-13 '""-13 

(3.8 



3.4.4 Canonical form of the equation for 7 

The only difference between equations (13.221) and (13. 29ft for uj and 7 is that, in general the co- 
efficient Tn 7^ 0. We can nevertheless perform the same change of variables as for u>, obtaining: 

Lemma 3.12. There exist coefficients dij(u>), < i,j < 3, and vector functions d^Ax^uo) such 
that the new function 71 defined as 

71 = 7 + d 20 z 2 + d 02 z 2 + d 30 z 3 + d 12 z 2 z + d 12 zz 2 + d 023 z 3 + z(f, d' 1Q ) + z(f, d' 01 ), (3.82) 

with dij = dji, is a solution of the differential equation 

7i = rn(a;)^ + ffl. (3.83) 

Furthermore T R satisfies the same estimate Ii3.28\) as T R . 

3.4.5 Summary of the equations in canonical form 

We summarize the main formulas of §3.4.H - §3.4.3I First we recall that 

f = g + h, g = P d T f, h = P c T f, h = k + k 1 + h 1 

where k and hi are defined in (13. 481) - (13. 491) . The equations satisfied by h and hi are, respectively, 
(see ( 13381 and (I33u1> ) 

h = C M h - P c T jE 2 [w, w] + H R , (3.84) 

and 

JH = C M h + H R . (3.85) 

The remainder H R is estimated in (13.441) and H R = H R + H', where H' is estimated in (13.551) 

and (EL56D . 

The second equation, given in Lemma F3.81 determines the evolution of u\: 

wi = Cl R , (3.86) 

where fi R is estimated in !3.59[ and uj\ and uj are related by (see (I3.57P ) 

uJi-uj = K(uj)\z\(\z\ + \\f\\ L ~ j ). (3.87) 

The third equation describes the evolution of Z\ (see (13.741) ): 

Z\ = i[izi + iKr\zi\ 2 zi + Z R , (3.88) 
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with the estimate (13.791) for Z R . ^From (I3.63p . z and Z\ are related by that 

Zl - z = K(uj)\z\ 2 . (3.89) 

The fourth equation is for the evolution of y = \zi\ 2 : 

y = 2Re (iK T )y 2 + Y R , (3.90) 

where ^ 

IXrI<NI^I- (3.9i) 

The negativity of Re iK? is a key point in the analysis and was proved in Lemma [3791 It is clear 
from the last equation that this condition gives a nonlinear damping effect in the evolution of 
the amplitude of the discrete mode - this is the crucial dynamical consequence of the Fermi 
Golden Rule, with obvious relevance to the large time behaviour of the solutions. 

3.5 A bound for \ut — uj\ and the initial conditions 

We will need a uniform bound for \ut — w(£)| on the interval [0, T]. This can be obtained by 
comparison with the function ui(t) as follows: 

\u>t — ^(t)] < | Wit — Wi(t)| + \ui\t — u>t\ + \w\{P) — ^Wl (3.92) 



: / |wi(r)|dr + ^(wr, \z T \ + ||/t||l~ ){\z t \ + II/tIIl^) 2 + U(u, \z\ + ||/|U» )(|*| + \\j u <_ 



,2 



< mBxK(u, \z\ + ll/IU- ) [J (\z\ 2 + Wfh^Ydr + (\z T \ + \\fa\\ L ~ g Y + (\z\ + ||/|| L - ) 

t 

where zt = z(T), fa = f(T). Using \u\ < \u \ + \u — ujt\ + \oj — ojt\, we have 
max 1Z(uj, \z\ + II /"Hl 00 - ) = 7^( m & x la; — out\, max (\z\ + II f \\l°°,_ ))■ 

We denote such quantities by the symbol ^(w, \z\ + ||/||l^ )• Then (j3.92p becomes 

T 

t 
As in fll,18p . we suppose the smallness condition: 

\z(0)\<eV*, WmWii <<**'*, (3-94) 

where e > is sufficiently small. Equation (13.891) implies \zi\ 2 < \z\ 2 + lZ(u, z)\z\ 3 . Therefore 

y(0) = \ Zl (0)\ 2 < e + K{u, \z \)e 3/2 . (3.95) 
^From the formula h = Pj>f — f + (P d — Pf 1 )/; we see that 

\\h(0)\\ Ll <ce s / 2 + 11^)1^-00^(0)11^. (3.96) 
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3.6 A bound for k\ 

Lemma 3.13. The function k\ defined in $-4-9\) satisfies the following bound: 

ii*ilu-„ ^ c I^°)I 2 (TtW - C T^W (3 - 97) 

Proof. Equalities (BUOJ) and (ESj) imply 

[C M (T)dT C T t+i f/3(T)dr(n+-n-) 

h = -e° fc(0) = -e o T; fc(0) (3.98) 

i 

Denoting i/ = f f3(r)dr, we obtain by expanding the exponential and using the idempotency of 

o 
projections (the Euler trick): 

e ^n± = u ± e iu + h t + pd 

Therefore 

e iKn+-n-) = ( n + e ^ + n T + p^)(n^ e -^ + n+ + P£) = n+e iu + Itpe -6 ' + P*. 

Note that Cy commutes with P^, hence 

t 

JC M (r)dr = e C T t( e iv u + + g-^jj- + p^_ (3 gg) 

Since /5 is a real function, both exponentials are bounded. Further, by (13.481) we have k(0) = 
a 2 oZ 2 (0) + a u z(0)z(0)+a 2~z 2 (0) with a {j defined in ( 13331) . (13331) . Therefore, the bounds (12151) . 
(EIED, (12X71) and assumption (13311) imply (13371) . D 

4 Large time asymptotics 

In this section we will make use of the dispersive estimates given in §2.2l to prove the asymptotic 
representation for the solution of (11.21) with initial data as in theorem 11.31 The idea is to fix 
an interval [0, T] and carry out the frozen spectral decomposition relative to the operator 
Ct = C(u>t) at time T, as described in §3.2L We then obtain bounds for certain majorants on 
this interval which are uniform in T, and thus make it possible to obtain the asymptotics for 
the solution in the limit T — > +oo. We will use the 1Z notation explained prior to (13. 7p and 
(I3.93P to express estimates and bounds concisely. 

4.1 Definition of majorants 

We define the quantities 

Mo(T) = max \uo T - u\ I — ^— | (4.1) 

V ' 0<t<T l l \l +£ tJ K ' 

Mi(T) = max \z(t)\{ —^— | (4.2) 

M 2 (T) = max \hAv~ [ — — | (4.3) 

21 



which will be referee! to in the following as "majorants", and denote M the 3-dimensional vector 
(Mo,Mi,M 2 ). Observe that (15771) . (15751]) and (137971) imply that / can be bounded in terms of 
these: 



L » = ^i(cj)(M 2 + ^M 2 ) 



so that control of M will allow complete control of the asymptotic behaviour of the solution. 
The goal of this section is to prove that if e is sufficiently small, M is bounded uniformly in T. 
This is done by first bounding the initial data, and inhomogeneous terms, in the equations in 
section [3.4.51 in terms of the ML,-, and then using the estimates for the homogeneous evolutions 
to self-consistently bound the ML,- in terms of themselves, uniformly in T > and e< 1. 



4.2 Estimate of the remainders and initial data 

Lemma 4.1. The remainder Yr defined in $3. 81\) satisfies the estimate 



\Y R \= TZ(e 1/2 M)- 



r 5/2 



:i 



(4.5) 



{l + et) 2 y/ei 

Proof. Using again the equality f — g + h — g + k + k± + hi, lemma 13.971 and the definitions of 
the ML,, the remainder Yr is bounded as follows: 

Yr = 1Z 2 (uj, |z| + ||/||l~,)|2| (M 2 + ||&i|U~ /3 + II h ih™ p ) 2 + \z\\ut - uj\(\z\ 2 + ||/Cl||L^ + ||/il||L^) 



+ki(ii&iiu^+iMU> 



K{e 



1/2 



e \V2 



£ \3/2 



l+et, 
e \ i/2. 



UqIvIi 



+ 



l + et 
e 



6 M2 , E 



l + et x (1 + t) 3 / 2 Vl + et 



e \ 3/2 \ 2 



e \3/ 2 
^— ) M 3 



l + et x (1 + t) 3 / 2 Vl + et 



l+et 



1 + t) 3 / 2 Vl + et 



3/2 



K(e 



l/2i 



r-5/2 



yr- ~ ^-^ - (i + im|) 5 , 



establishing (14.11) . 

Let us turn now to the remainder Hr = Hr + H' in equation (13.851) for h\ 
Lemma 4.2. The first summand Hr satisfies 

e \3/2 



D 



\H 



R\\Mp 



Tl{e 



1/21 



(l + Mif + e^l + IMI) 4 ). 



Proof. It follows from (I3.44p 



l# 



fl A< f 



ft 2 0+ bl + ||/||l 



l + et 



kl 3 + (M + I^t - ^l)(M 2 + 1 1 &i 1 1 1,^ + IIM-^ 



(4.6) 



+(M 2 + ll^ilU- + IMU~J 2 = 7e(e 



l/2i 



rat2 i 



l + et 

e \ 3 / 2 . „ \ ( e 



3/2 



.l + et x (1 + t) 3 / 2 Vl + et 
which implies (14. 6ft . 



l + et 
e 



1/2 



l+et 



l + et * (1+t) 3 / 2 Vl+et 



e \ 3/2 x 2 

M 2 



D 
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The second summand H' is represented as in (I3.55P where the A m are estimated in (j3.56p . 
For the A m we now obtain: 

Lemma 4.3. 

llAJlM^Kie^^-^y^^Ml + e^il + lMlf). (4.7) 

Proof. Estimate (13.561) implies 

HAnll^ = TZ 2 {uj,\z\ + \\f\\ L o. )\z\[\z\\u; T - u\ + (|*| + ||fci||i,« + ||/ii||i,~ ) 2 J 



e \V 2 



e \ 3 / 2 // e \ V 2 e ( e \ 3 / 2 \ 2 



ft(e 1/2 M)(^^) Mi 

\l + etJ l Wl + etJ (1 + t) 3 / 2 Vl+et 

which implies (14. 7p . D 

Now we estimate the initial data. Referring to the formulas at the end of §3.51 we have 

y(0) <e + K(e 1/2 M)e 3/2 = e(l + K(s 1/2 M)s 1/2 ) (4.8) 

||M0)||^ < re 3/2 + fti(w)|fa^-w|||/(0)|| £ « < ce 3/2 + n(e 1/2 M)e 2 M (l+Ml + e 1/2 M 2 ). (4.9) 

4.3 Integral inequalities and decay in time 

This section is devoted to a study of the system: 

y = 2Re(tK T )y 2 + Y(t), (4.10) 

h 1 = C M hi + H(x,t), (4.11) 

under some assumptions on the initial data, and on the inhomogeneous (or source) terms Y 
and H. Equation (14.101) for y is of Ricatti type, and is similar to (13.901) . while (14.111) is similar 
to (13.851) . First, for the initial data, we assume 

2/(0) < £2/o, WhWWMe < £ 3/2 h (4.12) 

with some constant yo and ho > 0. As for the source terms, we assume that 

|KW|£F (irSv3 (413) 

and that H(x, t) = Hi(x, t) + H 2 (x, t), where H 2 = ^2(Ct — 2ifj,Tm — 0)~ l A m , A m G X? with 

m 

the following bounds: 

— / e \ 3 / 2 

\\M\m, < An{Y^a) > ( 4 - 15 ) 

where the quantities Y, H 1 , A m are supposed to be given positive constants. All these assump- 
tions are motivated by the estimates of the remainders in §4.21 and by the final estimates we 
intend to prove on u, z, h and h\. Equation (I4.10p corresponds to equation (I3.90p and the 
assumption (14.131) on the source term has the form of estimate (14. 5p for the remainder Yr. Sim- 
ilarly, equation (14. lip corresponds to equation (13.851) and assumptions (I4.14l) - (l4.15p correspond 
to the inequalities (I4.6p ~ (14.71) . Finally, corresponding to (I3.75p . we work under the assumption 

Re iK T = -lrs\K T < 0. 
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Lemma 4.4. fl]Q Proposition 5.6]) The solution of ^-lOj ), with initial condition and source 
term satisfying \4-lty an d \4-13\) respectively, is bounded as follows for t > 0: 



\y(t) 



2/(0) 



l + 2lmK T y(0)r 
Let us consider equation (14.111) for h± 



— ( e \ 3 /2 
<cY[——\ , c = c(y ,lmK T ). 



(4.16) 



Lemma 4.5. The solution of ( |^.iJ| ), with initial condition and source term satisfying ft4.12\ ), 
Hl4-14\) an d \4-lty i i s bounded as follows: 



\\hi\\ L ~ <c{uj t ){^—^ (ho + Hi + J^A 

m 

Proof. The function hi(x,t) can be expressed as: 



JC AI (r)dr f fC M (r)dr 

hi = e° hi(0) + e° H(s)ds, 



(4.17) 



To establish (J4TT7I1 we use the representation ([3799]) and the bounds (127L4J) . fl2A5|) and (|27T7j> 
to deduce that 



MU°°, < ^^L IMP) II. 



t 



-P ~ (1 + t)3/2 ' 



c(cj t ) 



' (l + (t-s)) 3 / 2 




|-H"l(s)||^a + Pm(s)||A4a)^ 



< c{u T ) 



hr 



e \ 3 / 2 



1 + t 



ds 



i + (t"I S ))3/2(^(rr^) +E^d 



+ es 



3/2 



£ \3/2/ _ 



- c ^ T )(lTrf) ( h ° + Hl + 5Z Am ) ' 



since / *(1 + t - s)^ 3 / 2 (l + es)~ z / 2 ds < c(l + e£)~ 3 / 2 by [21 lemma 5.3]. 

4.4 Inequalities for the majorants 

In this section we estimate in turn the three majorants M ,M 1 ,M 2 . 
Lemma 4.6. The majorants M (T), Mi(T), andM 2 (T) safe/?/ 



,(T) = TZ(e 



1/2 T 



:i + 



£(1 + |M|) 2 



2 = ^(£ 1/2 M)fl + £ 1 / 2 (l + |M|) ; 



?e(e 



1/2 T 



;i + Mi) 3 + e 1/2 (l + |M|)^ 



a 



(4.18) 



(4.19) 



(4.20) 
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Proof. Step i) Using the equality f = g + h = g + k + k\ + hi and bound (I3.97P for k\ we have, 
using the notation defined prior to (13.931) 



M 2 + 



7^2 (w, \z\ + ||/||i»,)(|Fi||£» + \z\ + ||/h||l»j, 



3/2 



1/2* 



n{e 



n(e 1/2 M) 



e I e 

+ 



1+t) 3 / 2 \l + £t 



1+et 



1 + et 



so that 



1 + 



e 



I 2 + e^M-,) , 



\z\ 2 + \\f\\ L ^='R{e^WL J 



ij + e^l 



(4.21) 



Then (13T93D and flHD imply (|Qg|) . 

Sfep mJ Recall y = |^i| 2 satisfies ( I4.10p with Y = Yr, and Yr satisfies the inequality ( 14. 5 p which 
is exactly the condition 1JQ5J) with F = ft(e 1/2 M)(l + |M|) 5 . Using (KWf as well as (TCHD to 
bound the initial condition y(0), it follows that 



y < n{e 



1/2| 



5 \3/2 



1 + et 



Vl + et/ v 



Therefore 



M < j/ + ft(w)M <^( £ 



-1/2] 



1 + et 



e \ 3/2 

T — ) (l + IMD- 



1 + et/ 



from which (I4.19P follows. 

Step in) Let us now consider h\, the solution of (I3.85p . It has the form (14. lip with H = Hr = 
Hr + H', where Hr and H' identify respectively to Hx and H 2 . More precisely, using (J4.6P and 
( 14.71) . we have 

''l + M!) 3 + e 1/2 (H llurlv ^ 



Hx = TZ{e 1/2 U 

A m = TZ(e 1/2 M) (Mf + e l ' 2 {l + |M|) 3 J . 
Concerning the initial conditions, we know that /ii(0) = h(0). Thus by (14.91) 

h = c + ft(e 1/2 M)e 1/2 M (l + M 2 + e 1/2 M 2 ) 



Applying Lemma H~5| we deduce that 



M 



which implies ( I4.20p . 



^'Mrh) 



e \3/2 
+ 



l + M 1 Y + e y2 (l + \M\y 



D 
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4.5 Uniform bounds for majorants 

The aim of this section is to prove that if e is sufficiently small, all the Mj are bounded uniformly 
in T and e. 

Lemma 4.7. Fore sufficiently small, there exists a constant M independent ofT ande, such 
that, 

|M(T)|<M. (4.22) 

Proof. Combining the inequalities (14.18j) - (14.20j) for the Mj, one get a estimate of the form 

M 2 < n(e 1/2 M) [(1 + M x ) 8 + e 1/2 (l + \M\f 

Replacing M 2 in the right-hand by its bound (14.191) . we get an inequality in the form 

M 2 < n(e 1/2 M)(l + e 1/2 F(M)) 

where F(M) is an appropriate function . From this inequality it follows that M is bounded 
independent of e -C 1, since M(0) is small, and M(t) is a continuous function of t. □ 

Corollary 4.8. The function uj{t) has a limit uj + as t —> oo. Furthermore, the following 
estimates hold for all t > 0: 

\u+-u(t)\ < M^- t , (4.23) 

1/2 



K*>l £ M (TT7i) ■ (424) 

/ e \ 3 /2 
MU~, < M(— — J , (4.25) 



+ 
Kl + et, 

l -p ^ ciM) iT7f (426) 

Proof Since \ojt — oj{t)\ < M , then applying this result to \cu(ti) — o;(t 2 )|, we see that 

ui{t) is a Cauchy sequence. It thus has a limit, denoted u + and (14.231) holds. The next two 
results follow immediately, while the final one is a consequence of (14. 4p . □ 

4.6 Large time behaviour of the solution 

In this section we deduce from corollary 14.81 a theorem which describe a large time behaviour 
of the solution. Notice that in the decomposition / = g + h = g + h± + k + ki, a fixed time 
T has been chosen, and all the components depend on u(T). From the above proposition, we 
know that u(t) has a limit uo + as t — > oo. So we can reformulate the decomposition by choosing 
T = oo and ujt = uj + . Namely, let us denote P^ = P c (cj + ) and P^ = 1 — P^. We define 
/ = g + h where g = P^/ and h = P^/. We also decompose h = h\ + k + k\ where 

k = a 2 oZ 2 + On zz + flo2^ 2 , 



I: 

k 1 = -exp( C + (r)dT)k{0) 



where a^- = a,ij(u + , x) and C + = C(cu + ) + i (u(t) — u+ + j(t) ) ( 11+ — 11^ J . All the estimates 

previously obtained in §3.4.21 -^41 for finite T can be extended to T = oo and lot = ^+ without 
modification. Thus we have proved the following result: 
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Theorem 4.9. Let the conditions of theorem 11.31 hold. Then, for e sufficiently small, there 
exist C 1 functions uj(t), j(t), z(t) as in lemma l3J\ and constants w + eR and M > 0, such that 
for allt > 0: 

t 

j(/w(s)ds+7(t)) / \ 

i/j(x,t)=e° (ty iU {x)+z{t)u{x,u})+l{t)u*{x,u}) + f{x,t)), (4.27) 



and 

\u(t) - u + \ < M T -^ ) \z(t)\ < M( r ^) V2 , ||/(t)|| L . < M^— t) (4.28) 

so that uj + = lim u{t) E M. ^4 corresponding statement also holds for t — > — oo. 

5 Scattering asymptotics 

We have now obtained the representation ( 14.271) of the solution ip(x,t). In order to prove 
statement of theorem 11.31 it remains to: 

• construct asymptotic expressions for u(t), z(t), 7(t), which is done in section [57T1 following 
[2J, and then 

• to prove the existence of \l/±, and hence obtain the scattering asymptotics ( 11.191) ; this 
second stage amounts to the construction of the wave operator, and is carried out in §5.21 
by the study of some oscillatory integrals. 

5.1 Large time behavior of z(t), u(t) and j(t) 
We start with equation (j3.74p for zi, rewritten as 

z\ = ifjjZi + iK + \zi\ 2 zi + Z R 

with K + = K(u + ); by (13.791) the inhomogeneous term Z R satisfies the estimate 

Z R = n(e 1/2 M) f (1 + M 4 ) = 0{t~% t - oo. 

(1 + £t) 6 ' z y/et 

On the other hand, we have, from (14.8D and ( 14.161) . 

Given the estimate (14.241) for \z\, and obviously the same one for \z\\, we have 

i ' = ^' + '^ 1 + 2I.nL(0/ ' + Z - ^ = 0(0, «-~. (5.1) 

Since by assumption in theorem 11.31 y(0) = 0(e) we can write y(0) = eyo with y = 0(1). Let 
us denote 21m K + y = ek +J 5 = ±^ K + so that eK + y = §e& + (l — i5). The solution z\ of (15.11) 
is written in the form 



i f (i(ti)dti t a i f fi(ti)dti 

e ° 



z\ 



: , lfl lS) z 1 (0)+ e i (l + ek + s)^-^Z 1 ( S )ds 

1 + ek + t)^ ' L J 



(1 + 6M)^ M) 



+ ^R 
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where 

OO a 

/— i f u(t\)dt\ i ,. . „ 

e o {I + ek + s)^ {l ~ l5) Z^ds 

o 

and 

OO t 

f if»(t 1 )dt 1 /l + ek + s\^-iS) r7 . . , 
z R = - e° ___ Zi(s)ds. 



1 + e/c+t 
t 

Here /x(£i) = /z(a>(£i)). From the bound (15.11) on Z\ it follows that Zr = 0(t~ l ). Therefore 
z{t) = z x {t) + Oit- 1 ) satisfies 

t 

if n(ti)dti 

e ° 

z(t) = z + - — + 0(f~ ), t-»oo, 2+ = 2;oo(^+)- (5.2) 

(1 + ek + t)'^ ' 

^From these formulas for z(t), the asymptotic behavior of uj{t) and 7(£) can be deduced as in 
[21 Sections 6.1 and 6.2], leading to the following: 

Lemma 5.1. In the situation of theorem \4-£\ the functions u(t) and ^(t) have the following 
asymptotic behavior as t — » +oo: 



u{t) = lj + + --±±— + - ±—cos(2iL + t + b 1 log(l + ek + t)+b 2 ) + 0{r 3 / 2 ), (5.3) 

1 + ek + t 1 + ek+t 



^— + ±— cos(2//+i + 6i log(l + ek + t) + b 2 ) + 0(r 3 / 2 

f e£; + £ 1 + ek + t 

7 (t) 7+ + c+ log(l + efc+t) + ©(r 1 ), (5.4) 

where u + ,k + are as defined above, /i + = fJ,(ou + ), and q+,b + , b\,b 2 ,c + are constants. 

5.2 Soliton asymptotics 

Here we prove the statement (11.191) in our main theorem 11.31 To achieve this we look for the 
solution ip(x,t) to (II. ip . in the corresponding complex form 

^ = s + v + /, (5.5) 

where 

is the accompanying soliton, and 

v(x, t) = v(x, t)e ie(t \ v(x, t) = (z(t) + z(t))u!(x, w(t)) + i(z(t) - z(t))u 2 (x, u(t)). (5.6) 

We aim now to prove the complex form ( 11.20J) of the scattering asymptotics, by analysing /. 

Lemma 5.2. Ifif;(x,t) is a solution of /il.l\) . and f is as in ( 15.51) - ( I5.6f) . then 

if = -f" + R 

where 

R = 7(s + v) — iuod u {s + v) — i[(z — ifj,z)(ui + iu 2 ) + (z + ifiz)(ui — iu 2 )]e l 

S(x)[af + be ie Re (e~ ie f) + 0{\f + v| 2 )]. (5.7) 
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Proof. Multiply (11.101) by e l9 , we obtain — us = —s" — S(x)F(s) that implies 

is = —us — 7 s + iud^s = —s" — ■ys + iud^s — S(x)F(s), (5.8) 

By the equation Cm = ifiu, we obtain for the components u± and u 2 of vector u, 

j —u'2 + uu 2 — S(x)au 2 = i\xu\ 

\ —u'{ + uu\ — 5{x)[a + b]ui = —i\iu 2 . 

Therefore — v" + uv — S(x)[av + bRev] = —fi(z — z)u\ — ifi(z + ~z)u 2 and then 

iv = -( w + y)ve m + iwd u ve m + (i(z + i)u x - (i - i)w 2 )e ie (5.9) 

= — v" — 7V + iud^v — <J(x)[au + 6Rew]e je + i[(i — i/j,z)(ui + ra 2 ) + (z + i/Xz)(ui — iu 2 )}e ld . 

^From (II. ip . f)5.8p . (I5.9P we obtain for the remainder f(x,t) = ip(x,t) — s(x,t) — v(x,t) 

if = —f"(x,t) +7(5 + v) - iud w {s + v) - i[(z - ifiz)(u{ + iu 2 ) + (z + ifiZ)(ui - m 2 )]e 4e 
- 8{x)[af + 6e* e Re (e""/) +0(\f + v| 2 )] = -/" + R, (5.10) 

with i? as in (15771 D 

The function f(x,t) which is a solution of 15. 101 can be expressed formally as 

t 
fit) = W(t)f(0) + jw(t-T)R(T)dr 



00 00 

= W(t)(f(0)+ f W{-r)R(r)dr) - j W{t-r)R{r)dr (5.11) 

t 

= W(t)(f> + + r + (t), (5.12) 

where W(t) is the dynamical group of the free Schrodinger equation. To establish the asymp- 
totic behavior (jl.20p . it suffices to prove that 

0+ G C b (R) n L 2 (R), and ||r + (t)|| Ci(R ) nLa(R) = 0(t~ v ), t -► 00. (5.13) 

These assertions follow from the definition (15. 7p of the function R, and the following two 
lemmas. The first lemma studies the contribution to 4>+(x) and r + (x, t) from the terms in (15.71) 
involving 8(x) which is 0(t~ l ) as t — > 00 by (14.261) . 

Lemma 5.3. Let Tl(t) be a continuous bounded function of t > 0, with \H(t)\ < L and 
\tU(t)\ < Li. Then 

//* — ix / (4t ) 
W(-T)[6(-)U(r)]dr = / 6 / U(r)dr G C b (R) n L 2 (M) (5.14) 



and for v G (0, ~) there exists C = C(u, L , Li) > such that the function 

00 00 

//* ix /4(£ — t) 

W(t-r)[5(-)n(r)]dr= / 7 =n(r)rfr (5.15) 

7 v / 47ri(t - r) 

satisfies \\r(;t)\\ Cb (B)ni?(8) < C{l + t)~ u . 
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Proof. The CVproperties follow from formulas (j5.14p and (15.151) (in fact with v = 1/2). To 
prove the L 2 -properties, let us change the variable to t = 1/u to get: 



oo 



i ' lux2/i r,(u) du = -^=^ u ^ x yMu)), V(u) = Ii{l/u)/u^\ (5.16) 

where J r u ^,^(f(u)) = /(£) indicates the Fourier transform with argument £. By the assumptions 
on II we have \f](u)\ < Lqu~ 3 ^ 2 as u — > oo, and 1 77(11) | < Liu" 1 ^ 2 as u — » 0. Therefore 
77(14) G L P (M) for 1 < p < 2. It follows from the Hausdorff- Young inequality for the Fourier 
transform that (f> G L 9 (R) for q > 2 as a function of y = x 2 , i.e. 

00 



and hence <fi G L 2 (M), since it is already known to be bounded and continuous. It remains 
to prove the decay rate of r(x,t) in the norm L 2 (IR). Let us represent the function r as 
r(x, t) = W(t)p(x, t), where 

00 1/t 

p(x,t) = fw(-r)[8(-)U(r)]dr = -=L= /V^ 4 rj(u) du = -^F Vr ^, A {C t {u)r,{u)). 

J \/-4:iTi J V-2^ 

t 

Here Ct( M ) is the characteristic function of the interval (0, 1/t). As above p is bounded, but also 
since 

l/t l/t 

\\(t(u)r)(u)\\ LP = ( / |r/(n)| p rfn) 1/P < L x ( f U - p/2 du) ^ < Cr^ , 1 < p < 2, t > 1. 



the Hausdorff- Young inequality implies that for any q > 2, and in fact as t — ► 00: 

lb(x,t)(i + |x|) 1 /«||L,<c(r 1 ^ a ) ) 

for some constant C = C(Li,p), for g _1 + P" 1 = 1. The Young inequality then implies that 

\\ P (z,t)\\* < IIpOM)(i + M) 1/9 IMI(i + Wr 1/9 lk = O(r^), i = - + -, 

2 q r 
if r > q. To have r > q, we must take g < 4, or equivalently p > 4/3. Hence, we have 

„ = i^/?< 1/4 . □ 

p 

The second lemma studies the contribution to <fi+{x) and r + (x, t) from terms without 5(x) 
in (15. 7ft . Consider the expansions ( 13.22(1 . ( 13.29(1 . ( 13.33(1 . for u)(t), j(t), and i(t) — ipz{t): the 
main (quadratic) parts of these contain the terms £+(£), £+(£), 2 +(£)z+(£), which are 0(t~ l ) as 
£ — > 00. The remainders are 0(£~ 3 / 2 ), and it is straightforward (from the unitarity of W) to 
bound the contribution of these to 0+ in Cj n L 2 , and to check that these contribute 0(t~ l l 2 ) 
to r + in C*6 fl L 2 . Thus, without loss of generality, we may replace u(t), 7(t), and z(t) — ipz(t) 
by the main quadratic parts. We first show how to treat these terms with the phase 6(t) in 
(15.71) replaced by <£>+(£) = w+t, and then consider the general case in lemma 1531 
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Lemma 5.4. Let U(t) be one of the functions z 2 + {t)e l ' p+ ^ , \z + {t)\ 2 e lip+ ^ or z 2 l _(t)e' l ' p+ ^ , where 



gi/i+t 



z + {t) = — — —- j-, with ek + > 0, tp+(t) = u>+t, and let ifj(x) be a bounded continuous 

function of x 6l, such that x 2 ijj(x) G L 2 (IR). Then 

oo 

/ n(r)W(-T)ipdT G C b (R) n L 2 (R) (5.17) 

o 

and for each v G (0, 3 ~ 4 ) there exists C v > such that 

n(r)W(t - r)Hr\\ Cb{R)nL2{R) < CjT\ (5.18) 



Proof. Since ||W^(£)^||c,, = 0(t l ' 2 ) then C b - properties are evident, in fact with v = 1/2. It 
remains to prove the L 2 -properties. Since ip £ -^ 1 (M) H L 2 (IR), we have 

W(t)V = ^L=J e^-y\ 2 l u ^{y)dy 

ix^ /At p ix^ /At r 

~"*y/^(y)dy + e -—= \ e -i*y/2t( e iy 2 /*t _ l)^(y)dy 



y/47rit J ^/Airit 

Ax 2 1 At „ix 2 /4t 

= ^ x /2t) + —=Mx/2t), (5.19) 

\J lit V lit 

where ipt{y) — { e%y ^ U ~ l)i*{y)- For t > 1 we have, using \e %e — 1| < 9, 

-±=\\M-/it)\\ L . = \\U)h> = \\M-)\\v = (J\(e iy2/u -my)\ 2 dy) 1/2 < ^||i/V(i/)I|l» < f 

Therefore, using the fact that |n(r)| < (1 + ek + r)~ l , we deduce that 

oo oo 

— ix J At /* — ix I At 

n(r) e y , j) T {-x/lr)dT G L 2 (R), and / n(r)^^^ T (-x/2r)rfr = 0(t _1 ) 
V-2zr 7 V-2rr 

o t 

in the norm of L 2 (M). Hence, to prove (15.171) it suffices to check that 

oo 

— 23; At 

n(r) 6 . . ^{-x/lr)dT G L 2 (M), (5.20) 

V-2zr 
o 

and to prove (I5.18P it suffices to check that 

00 ■ 2 I A 

— ix J At 

Y[{r) e - T ==^{-x/lT)dT = 0(t~ v ), (5-21) 

t 

for appropriate v > 0, in the norm of L 2 (IR). 
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First, let us prove fl5T20l) and fl5T2TD for IT(t) = z 2 At)e i(p+ ® . Note that \d T ^(x/2r)\ < C 



T 2 



because xip(x) G L 1 (M), and then ijj' is bounded function. Therefore we obtain (I5.20p with help 
of integration by parts: 

\<j>(x)\ = iT-^ ^(~x/2t) AT 2 / -^/4r-M( W+ +2 W .)rN d I (533) 



oo 



2r(l + e£; + r)(x 2 + 4t 2 (cj+ + 2/j+)) 
^(-x/2t)r 3 / 2 



2/ 



- C X ! r( (l + efc + r)(x 2 +4r 2 ( W+ + 2/i + ))' )| " T+~R ' 

Similarly, (I5T2T1) follows with v = 1/2. For the function U(t) = \z + (t)\ 2 e^ +{t) the proof of <K2Hh 
and (j5.2ip is similar. 

Second, we consider the case IT(£) = ~z 2 ,(t)e l(p+ ( t > which is more difficult, because the factor 
x 2 + 4r 2 (u; + — 2/i + ) vanishes for r = t* := \x\/2^2fi + — u+. We consider only positive values 
of x since the negative values can be considered similarly. Let us choose 1/2 < p < 1 and 
1 < q < Ap — 1. Then for large x we have 

0<t* -x p <t* <t* +x p <x q . 

Let us represent (0, oo) = J\ U J2 U J3, where Ji = (x q , 00), J2 = [£* — x p ,£* + x p ], and 
J3 = (0, 00) \ J\ U J2. Then = 0i + 2 + 03, where 0$ is obtained by integrating the integrand 
in (15.201) over r G Jj, i = 1, 2, 3. For 0i and 02, it is immediate, without the need to integrate 
by parts, that 

00 t*+xP 

|0i(x)| <cj^ = cx^/ 2 , |0 2 (x)| < ^ J dr< cx-W 2 -?\ 

xl t*-xP 

Therefore 0i, 02 G £ 2 (1, 00) since 3/2 — p > 1/2. Next observe that 

\x — 2ta/2/i + — lu + \ > cx p and r < x q for all r G J3. (5.23) 

Therefore for large x we use integration by parts to obtain 

f ( CT 3 / 2 ~ \ 

\M*)\ < ™' + ^ iar( (1 ; tk+T)[x + 2rV2 ^, + _ , +)( , _ 2t ^ + _ g^ V-(-x/2r)j| dr 

s CT "' + c /(rT$^ + c /(TT7^ £ ^ + CT8/ ^ (5 - 24) 

J3 J3 

(The boundary terms arising in this integration by parts can be estimated to be < c r 1 ' 2 x~ p /(x+ 
t) < ct- l l 2 x- p .) Since 2p - q/2 > 1/2 it follows that 3 G L 2 (l, 00). 

Similarly, we can estimate p. Again we assume, without loss of generality, x to be positive. 
For bounded x, say < x < I, I > we can just estimate the integral (15.211) directly as 

, M , f^ CdT 

p(x,t) < r \ <<*-*, x<l. 

Jt vn 1 + efc + r ) 

The decay of \\p(-, £)||l 2 follows from this together with a bound of the form \p(x, t)\ < ct~ u \x\~ a , 
with a > -z, which is valid for large t and \x\ > I 3> 1. To obtain such a bound we decompose 
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the integral (15.211) as p = Y^i=i Pii where Pi(x, t) is obtained by restricting the integral in (15.211) 
to r G Ju = Ji fl [t, oo). 

To bound p\ there are two cases. Firstly, if x q > t then |pi| < c J°^ r~ 3 ^ 2 dr < cx~ q l 2 . 

Secondly, if x q < t then |pi| < c J t °° T~ 3 l 2 dr < ct' 1 ^ 2 . In both cases \p\\ < ct~ v \x\~^~ u " >q . 
To bound P2, we notice that r > ex on J2t, and also r > t, and then just estimate 

|p 2 (M)| < c [ f T r ^ < cminr- 3/2 |x| p < cr^lxr 3 ^^. 



To bound p 3 , notice that for r G J-^t the inequalities (15.231) hold so that it is possible to 
integrate by parts (as in (j5.22p above) since the denominator which appears is bounded below. 
As with 03, the boundary terms arising in this integration by parts are < cr 1 ^ 2 x~ p /(x + r) < 
ct~ 1 l 2 x~ p . Also, as in ( 15.24ft . the integral which remains after this integration by parts can be 
bounded as < ct~ 1/2 x~ p + cx q/2 ~ 2p < ct~ 1/2 x~ p + ct-"x q/2 ~ 2p+uq , since t < r < x q in J 3t , so that 
we may assume t < x q in the estimation of p%. Therefore, in conclusion we have the following 
estimate for x > I ^> 1 and t large: 

\p(x,t)\ < Ct-"(x- p + X- 3/2+P+U + x -2p+(l/2+«/)9 + x -(l/2-u)g^ ^^ 

This shows that p(x, t) is square integrable and (15.211) holds if 

p>l/2, 3/2-p-u>l/2, 2p-(l/2 + u)q>l/2, (1/2 - v)q > 1/2 

The conditions on the exponents can be written equivalently as 

1 4o- 1 

<q<— , l/2<p<l-u. (5.26) 



l-2z/ l + 2z/ 



Such g and p exist if < v < (3 — \/5)/4, since then |-|^ > y^; > so choosing p = 1 — i 

for e small and positive, and g in the interval (15.261) will work. □ 



c>\ 



Lemma 5.5. The conclusions of lemma 5^4 o,re valid if <p + zs replaced by 9 satisfying 9{t) — 
<P+(t) = (^(t -1 ) /or large r. 

Proof. The estimates in the proof of lemma 15.41 which involve estimating the integral of the 
absolute value are completely unaffected by change of phase, so it is only necessary to re-assess 
the argument involving integration by parts, i.e. the treatment of 03 and p%. For example, in 
the more difficult case when IT(£) = ~z 2 + (t)e %e ( t \ we proceed as follows in the estimation of P3. 
Write = 6 — <f> + , and integrate by parts exactly as before, leaving along the e 1 ^ factor: this 
factor then carries throught to the integrand after the integration by parts, and we need to 
bound: 

M ' T \(l + ek + T)(x + 2T^/2p-^=ZJl)(x-2T^2 l i + -u J+ )^~ X ' 2T >r dl ~ 

(The treatment of the boundary terms is unaffected since \e l ^\ = 1.) But since by assumption 

0(t) = 0(t~ 1 ) the extra contribution to the integrand can clearly be estimated for large x, r 
in the same way as the term arising from differentiation of r 3//2 , and so f 1 5 . 2 5 1) still holds, as 
required to complete the proof. □ 

Remark 5.6. The t — > — oo case is handled in an identical way. 
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A The eigenfunctions of the discrete spectrum 

Here we find the function u = u(u) satisfying Cm = Am, where A = ifi. Using the definition of 
the operator C, we obtain 



/ 



-A 



d 2 



dx' 



+ U) 



d 2 
\dx~i-" 
If x 7^ 0, the equation flA. II) takes the form 

d 2 



u = S(x) 



A 



a 
-a-b 



u 



(A. i; 



dx 2 



-A 



U) 



dx' 



+ UJ 



-A 



\ 



u = o, x^O. 



(A. 2) 



General solution is a linear combination of exponential solutions of type e lkx v. Substituting to 

COD , we get 

-\ k 2 + u\ 
,2 , , x )v = 0. (A. 3) 



-Jfe a - c^ 



-A 



For nonzero vectors w, the determinant of the matrix vanishes: A 2 + (k 2 + u) 2 = 0. Then 
k\ + u; = =pA. Finally, we obtain four roots ±fc±(A) with 



fc-i-(A) = V— lj =F *A, 



(A. 4) 



where the square root is defined as an analytic continuation from a neighborhood of the zero 
point A = taking the positive value of Im^—u at A = 0. We choose the cuts C + in the 
complex plane A from the branching points to infinity. Then ImA;±(A) > for A G C \ C±. 
It remains to derive the vector v = (v\, V2) which is solution to ( 1A. 31) : 



V 2 



kl +uj ±i\ 



A 



A 



-vi 



v\ = ±iv\. 



Therefore, we have two corresponding vectors v± 
dent exponential solutions. 

' 1 



±i 



and we get four linearly indepen- 



v+e 



±/A:+.l' 



,±(A'+.f 



f_e 



±ik—x 



^zizik—x 



Now we find the solution of (1A. II) in the form 

u = Ae ik+lxl v+ + Be lk ~ lxl v_. 
At the point x = we have a jump: 

' a + b 



«-w-('r:) 



u(-0) 



Substituting (I A. 51) . we get 



2ik + Av + + 2ik^Bv^ = -M(Av + + BvJ), M 



a + b 
a 



(A. 5) 



(A. 6) 



(A. 7) 
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Note that 

Mv + = av + + f3v- , , 

, where a — a + -, /3 = -. 

Mv- = av- + (3v + 

Then (I A. 7p becomes 

(2zfc+ + a) A + (3B = 

/M + (2ifc_ + a)5 = 

The determinant D = (2ik + + a)(2ik^ + a) — /3 2 vanishes for A = i\x since i\x is spectral point. 
Therefore, set A = 1, and obtain 

u = e ik+lxl v+ -^ e ifc - |a V (A. 8) 

+ 2ik„ + a v ; 

Note, that 2ik- la^O. Indeed, if 2i/c_ + a = 0, then /3 = 0, a = a. and 2i/c_ + a = 
-2-y/u; + n + a = - sj a 2 + 4/x + a < 0. 

Since both fc + = a/— a; + // and fc_ = ^/— a; — ju are purely imaginary, the first component «i 
is real, while the second one «2 is imaginary. It is easy to prove that u* = (u\, —112) is the 
eigenfunction associated to A = —ifi. 

B The eigenfunctions of the continuous spectrum 

Consider A = iv with some v > u. 

I. First we find an even solution u = t + of equation flA. 1|) in the form 



r+ = (Ae lk+lxl + Be~ lk +W)v + + Ce ifc - |a V. (B. 1) 

At the point x = we have, similarly (I A. 6j) and (IA. 7j) . 



2ifc+(A - 5)*;+ + 2ik.Cv. = -M(A + 5)t; + + Cv. 
which equivalent to the system 

2tk + (A -B) = ~a(A + B)-f3C 

2ik^C = -/3(A + B)-aC 



(B.2) 



2ik 
If 2ik_ + a = 0, then D = -f3 2 , (3 ^ 0, A + B = 0, and C = -±{A - B). Put A = D we 

obtain B = -D, C = A(3ik + . 

If 2z/c_ + a^0, then from the second equation of (IB. 2D we get 

C = L—U + B). (B. 3) 

a + 2ik_ ' K ' 



Then, set A = D, we obtain from the first equation of (IB. 2j) 



, _ . -a(a + 2ik)+(3 2 / — 
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Solving this equation, we get B = —D and then C = 4[3ik + . Finally, we obtain 

T+ = (De ik+]x] - De~ ik+]x] )v + + A^ik + e ik -^v.. (B. 4) 

II. It is easily to check that an odd solution u = s of equation (1A. II) is 

s = ^(e ik+x - e~ ik+x )v + = sm(k+x)v + . (B. 5) 

For A = iv with v < —uo we have similarly 

r_ = (De ik ~ lxl - De~ ik -W)v- + Apik.e tk+lxl v + , s = sin{k-x)v- 



C Proof of Proposition 12.3 

First we prove the following lemma 

Lemma C.l. Let K(£) be the integral operator with the kernel 

K{t,x,y) = J KWV^KW'V) dU} (C.l) 

\u—vo\<5 

where f is a smooth function, f{v,x,y) = for \v — u \ > 5, and 

\d»f(v,x,y)\<c(N)(\x\ + \y\) N , N = 0,1,2,... (C. 2) 

Then for any a > N 

\\K(t)\\ Ba = 0(t- N ), t^oo, with a>N. (C. 3) 

Proof. For z = \x\ + \y\ > 5 we split the integral in the right hand site of (1U. ip into two part 
and obtain 

\ vl . \\ ^ f J( u ,x,y)- f(u ,x,y) f f(v,x,y)- f(v ,x,y) 
\K(t,x,y)\< I \ dv + / dv 

J V - Vq J V - V 

\u-v \<± ±<\v-v \<6 

5 

< -c{l)z + 2c(0) / — = c(l) + 2c(0)[logz + c]. 

1/z 

Hence, ( 1C. 31) follows with N = 0. Equality ( 1C. 31) for arbitrary iV > 1 can be obtained by 
similar way, using the integration by parts and inequalities ( 1C. 21) . D 

Proof of Proposition [Ol The operator e ct (C — 2i/i — O)^ 1 admits the Laplace representation 

e ct (C - 2ia - or 1 = / e xt R(X + 0) dX R(2iu + 0). 

2m J 
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Let us apply the Hilbert identity for the resolvent: 



R(X 1 )R(X 2 ) = [fl(Ai) - R(\ 2 )}, ReA fc > 0, k = 1, 2 

Al — A2 



to Ai = A + and A2 = 2«/x + 0. We obtain 



CDC 






r/A 



-JOO 



(OG 



W filA + 01-^+0) ^ / e , i(1 _ (A)) ^(A + 0)-fi(2^H-0) 

27ri y V ; \-2ifi, 2ni J v si ;; \-2iy, 

—ioo C+UC_ 

I e , (1 _ c(A)) J KA + o)-^ + o) rfA = Ki(() + KM + K3(() _ 



2-iri 

(-ioo,ioo)\(C+UC_) 

where C(A) G Cg°(iR), C(A) = 1 for |A - 2z/i| < 5/2 and C(A) = for |A - 2i(i\ > 5, with 
< 5 < 2fi — uj. By Lemma fC . 1 1 with N = 2, we obtain that 

\\Ki{t)\\ Bp =0{r 2 ), *-oo, 

since /? > 2. The bounds 1JTT2J) for /(z/) = i2(A + 0) follow from formulas 1J55J1 - (1231) . 

For the operator K2(i) we can apply the arguments from the proof of proposition 12.21 and obtain 

\\K 2 (t)\\ B$ = 0(t" 3/2 ), t^oo. 

Further, the integrand in K 3 (t) is an analytic function of A 7^ 0, ±i/i with the values in Bp for 
(3 > 0. At the points A = and A = ±«/x the integrand has the poles of finite order. Hoverever, 
all the Laurent coefficients vanish when applied to P c h e X c . Hence for K^i) we obtain, twice 
integrating by parts, 

||K 3 (t)P c fc|U~ <c(l + t)- 3 / 2 ||/i||^, 

completing the proof. □ 



D Proof of Lemma 13.5 



We use the following representation (see pQ): 
P C T = -^ / (R(A + 0) - R(A - 0)) dA 

C+UC_ 



— /(R(A + 0) - R(A - 0))dA + — /(R(A + 0) - R(A - 0))d\ = n+ + n~.(D. 1] 
27ri J 2iri J 



c+ c_ 

Let us decompose the resolvent given by (\2AL and (12.5J1 ) as 
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R(A, x, j/) = T(A, x, j/) + P(A, x,y) = ^2 A k(\ x, v) T k, (D. 2) 



fc=i 
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where 

ik+\x-y\ _ ik + (\x\ + \y\) in — 2k , „ , , ,s 03 

1 4fc+ ' 2£> 3 2D 

4, - _i£ p ifc- *L**+ » 4 K _ %a + ZK + r ik-(\x\+\v\) A _ _^ e _ IT) X\ 

Aa ~ 2D e 6 ' A5 ~ 2D 6 ' Aq ~ Ak„ AD - 6) 

and 

ri=r2= ") ,r3= -0' T4= (- : ,T5=T6= (- (D - 4) 

For the matrices Tk we obtain 

nj' 1 = irk, k = 1,2,4, Th]' 1 = -irk, k = 3,5,6. 
The terms A\ and Aq disappear when substituting (ID. 2[) - (lD. 41) into (ID. lj) since 



Ai(A + 0) - Ai(A - 0) = 0, Ag C_; A 6 (A + 0) - A 6 (A - 0) = 0, Xe C+, 
and we get 

P^- 1 - *( n T - U t) = ^- f[2(M* + 0) - 4»(A - 0))r 3 + 2(A 5 (A + 0) - A 5 (A - 0))r 5 ] rfA 
+— f[2(A 2 (X + 0) - A 2 (A - 0))r 2 + 2(A 4 (A + 0) - A 4 (X - 0))r 4 ] dA. 

27T2 J 

C- 

Let us consider only the integral over C + ; the integral over C_ can be dealt with by an identical 
argument. For A G C + we have: fc + = \/—uj — iX is real, and k + (X + 0) = —k + (X — 0) while 
fc_ = V - ^ + ^A is pure imaginary with ImA;_ > and A;_(A + 0) = A;_(A — 0). 

Note that ^(A, x, y) for A G C + exponentially decay if \x\, \y\ —> oo and smallest exponential 
rate of the decaying is equal to (2a;) 1 / 2 . 

It remains to consider the integral over C + with integrand A3 (A + 0) — A3 (A — 0). We change 
variable: £ = \/— u> — iX for the first summand and £ = —\/—uj — iX for the second. Then we 
get 



00 



I(z,y) = I (A 3 (X + 0)-A 3 (X-0)) dX = ^ e -V^\v\ { __ _ ) dX 



+00 



, y/2u,+C 2 \y\ p i<\x\ 



where 



D± = a 2 — (3 2 ± 2iaV— UJ — iX — 2a\/u) — iX =F 4zv— w — iXVuj — iX 
and 

D(C) = a 2 - (3 2 + 2ic< - 2a v / 2^ + C 2 - 4< v / 2^ + C 2 - 
Writing e"^' <i£ = ^(ie 4 '''^, and integrating by parts, we get that 



\I(x,y)\ < Ce- V2 ^ M (1 + \x\)~ m , Vm G N, 
completing the proof of lemma 13.51 
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E The Fermi Golden Rule 

In this section we study further the condition embodied in the Fermi Golden Rule ( 11.171) , 
showing that it holds generically in a certain sense: in particular, if a(-) is a polynomial function 
then, generically, the set of values of C for which ( 11.171) fails is isolated. To start with let us 
use the formulae in Appendices A and B to express the Fermi Golden Rule in terms of C and 
a(C 2 ) and its derivatives. By (IB. 41) 



T + (2ifj.) \ x=0 = (D - D)v + + 4(3ik + v„ = -4ik + (a + 2ik-)v + + A(3ik + v„ (E. 1) 

/ \ ( k + 1 \ .„.. a + 2ik_ , — — 

a{Kv + + V-) = a [ . f ^ ^ ], a = 4pik +J k = , k± = y / — uo ± 2p. 







Represent E 2 [u,u] = S(x)E 2 [u(0),u(0)], where 

E 2 [u(0), «(0)] = a'(C 2 )(u(0), u(0))*(0) + 2a"(C 2 )(^(0), u(0)) 2 ^(0) + 2a'(C 2 )(^(0), u(0))«(0) 

BydAH) 

p + 1 \ 2i/c_ + a , 

, fc_ = V— u - p. 



Therefore (u(0), w(0)) = (p + l) 2 - (p - l) 2 = 4p and 



E 2 [u(0),u(0)}=a'(C 2 )4:p 



° ^+2a"(C 2 )C 2 (p+V 2 ^ ° 







l) [C-\) 



Using (IE. lj) and (IE. 21) . we obtain 



(E.2) 



(r+(2i», £ 2 [w, u]) = <r(« + 1) a'4pC + 2a"C 3 (p + l) 2 + 2a'C(p + l) 2 + <r(« - l)2a'C(p - 1) 
Therefore, the Fermi Golden Rule fll. lTj) is equivalent to the condition 

;« + i)2p + (« + i)(p+ 1) 2 + (« - i)(p 2 - i)l + o"(« + i)c 2 (p + 1) 2 ^ o 



or 



with 



a" ^ 



2a'(2fi;p + 2p + fi;p 2 + l) 
C 2 («; + l)(l + p) 2 



(E.3) 



a + ~ - 2V^ + 2p 



6/2 



a + 1 - 2v^7T7^ a 2 b b 2 

6/2 ' 4 ' P 4 V 4 



i/a 



z/a 



Let v = — . Then u = — \J a? — v 2 a? = y/l — v 2 and 

2a 2 2 



a + ua- 2J\ + ua 2 VT 



v- 



l+i/ 



1 + 4i/Vl - ^ 2 



i/O 



a + z,a-2 v /^ + ^-v / T^^ 1 + „ _ y/i + 2^1 - v 2 



mi 
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are also functions of v only. Thus we conclude that there is a function c = civ) such that the 
Fermi Golden Rule holds if and only if 



C 2 ' a(C 2 ) ' 

Notice that the function c(V) is algebraic. It follows from this that if the function a(-) is 

polynomial, or even real analytic, then generically the Fermi Golden Rule holds except possibly 

at a discrete set of values of C. To see this observe that if the set of points where F{C 2 ) = 

c(u)a'(C 2 ) 
a"(C 2 ) — vanishes has an accumulation point, then F must be identically zero 

c(u)a'(C 2 ) 
since it is real analytic. But the condition a"(C 2 ) = — is a second order ordinary 

differential equation which determines the function a(C 2 ) given its value a(Cg) and that of its 
first derivative a'{Cl) at any point C = Cq. Clearly a generic polynomial function a(C 2 ) will not 
satisfy this equation, and so the set of points where the Fermi Golden Rule fails cannot have any 
accumulation points generically. This is also true for real analytic a in the following sense. Fix 

c(is)a'(C 2 ) 
any Cq, then there is a two parameter family of functions a(C 2 ) for which a"(C 2 ) = — = ; 

(this family of exceptional functions is parameterized by a(Co), a'(Co)). If a(C 2 ) is not one of 
these functions then the set of values for which the Fermi Golden Rule fails is at most a discrete 
set. 
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